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PREFACE 

The  results  reported  herein  were  developed  for  the  Arnold  Engineering 
Development  Center,  Air  Force  Systems  Command,  by  Nielsen  Engineering  &  Research, 
Inc.,  under  Contract  F4060O-77-C-0O08-P000ft|.  The  Air  Force  Project  Manager  for 
this  contract  was  E.  R.  Thompson,  AEDC/DOT.  This  report  covers  the  work  done 
during  the  period  September  1,  1978  to  December  31,  1979.  The  reproducibles  used 
in  the  reproduction  of  this  report  were  provided  by  the  author. 
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1.  INTRODUCTION 


The  flow  over  the  afterbody  of  jet  aircraft  is  often  char¬ 
acterized  by  separation  which  increases  the  drag  and  decreases 
engine  performance.  For  subsonic  flows  the  separation  occurs  in 
an  adverse  pressure  gradient  which  causes  the  flow  to  break  away 
from  the  wall  and  form  a  local  area  of  reversed  flow.  The  point 
along  the  wall  dividing  the  forward  and  reversed  flows  is 
defined  as  the  separation  point.  while  an  exact  flow  solution 
for  the  separated  region  can  theoretically  be  obtained  from 
finite  difference  calculations  of  the  fully  turbulent  compres¬ 
sible  Navier-Stokes  equations,  the  complexity  and  expense  of 
such  calculations  at  the  present  time  make  that  approach  un¬ 
realistic  for  engineering  requirements.  Consequently,  calcula- 
tive  techniques  have  been  developed  which  consist  of  patching 
together  approximate  models  of  the  flow  in  the  various  regions. 

All  of  the  methods  to  be  discussed  herein  apply  to  axisym- 
metric  bodies  at  zero  angle  of  attack.  In  the  method  of  Presz, 
et  al.  (ref.  1) ,  the  separated  region  is  represented  by  a  dis¬ 
criminating  streamline  which  determines  the  effective  shape  of 
a  solid  body  upon  which  an  attached  boundary  layer  is  calculated. 
Control  volume  analyses  locate  the  separation  point  and  deter¬ 
mine  the  shape  of  the  discriminating  streamline.  Boundary 
layers  and  exhaust  plume  entrainment  are  modeled  using  integral 
methods.  The  technique  is  applicable  to  subsonic  flows  with 
Mach  numbers  less  than  that  for  the  onset  of  shock  induced 
separation,  typically  a  free  stream  Mach  number  of  about  0.9. 
Another  approach,  developed  by  Cosner  (ref.  2) ,  employs  an 
integral  boundary  layer  method  for  the  entire  body  including  the 
separated  region.  Empirical  relations  are  used  to  extend  the 
calculations  across  the  separation  and  reattachment  points  for 
bodies  with  solid  plume  simulators.  Exhaust  plume  entrainment 
is  calculated  by  simply  extending  the  boundary  layer  analysis 
to  the  plume  boundary  treated  as  a  moving  wall.  Comparisons 
with  data  indicate  the  method  is  not  successful  for  bodies  with 
extensive  separation.  The  mast  sophisticated  approach  to  the 
problem  of  exhaust  plume  entrainment  is  that  of  Dash  et  al . 

(ref.  3)  wherein  finite  difference  methods  are  used  for  the  flow 
calculations  in  all  regions  of  the  flow.  Separated  flow  on  the 
nozzle  boattail  is  not  considered.  However,  studies  decribed  by 
Wilmoth  (ref.  4)  indicate  that  the  method  correctly  accounts  for 
exhaust  plume  entrainment.  The  present  author  developed  a 
method  (ref.  5)  which  combines  a  finite  difference  inviscid  flow 
method  with  integral  methods  for  the  boundary  layer  and  exhaust 
plume  mixing  layer.  The  approach  employed  a  displacement  sur¬ 
face  calculated  from  the  viscous  layer  which  forms  an  effective 
boundary  for  the  inviscid  flow  and  is  compatible  with  the 
resulting  inviscid  pressure  distribution.  For  separated  boundary 
layers,  the  displacement  surface  was  assumed  to  take  the  form  of 
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a  conic  section.  The  method  also  included  an  approximate  one¬ 
dimensional  exhaust  jet  model  including  entrainment  of  the 
boattail  boundary  layer.  The  method  was  found  to  provide  an 
accurate  engineering  prediction  method  for  boattail  flow  fields 
with  moderately  underexpanded  exhaust  flows  and  for  bodies  with 
solid  plume  simulators  for  free  stream  Mach  numbers  less  than 
that  for  which  shock  induced  separation  would  occur  on  the  boat- 
tail. 


In  this  report  improvements  and  extensions  of  the  method 
of  reference  5  are  discussed.  The  improvements  are  the  removal 
of  the  assumption  that  the  separated  displacement  surface  is 
conical  and  the  use  of  a  more  accurate  inviscid  jet  model  based 
on  the  method  of  characteristics.  The  calculation  method  is 
extended  to  include  shock  induced  separation. 


2.  SUBSONIC  INTERACTIONS 


The  particular  boundary  layer  and  inviscid  external  flow 
method  used  in  the  calculations  were  described  in  references  5 
and  6  and  will  not  be  described  in  detail  here.  The  boundary 
layer  method  is  an  integral  method  which  can  be  applied  in  a 
direct  mode  (boundary  layer  edge  velocity,  prescribed)  or  an 
inverse  mode  (boundary  layer  displacement  thickness,  6*  pre¬ 
scribed)  .  The  inviscid  flow  method  used  is  the  South- Jameson 
(ref.  7)  theory  as  implemented  by  Keller  and  South  (ref.  8)  with 
minor  modifications  to  accomodate  the  iterative  interaction 
procedure. 

In  the  method  described  in  reference  5,  the  first  step  in 
calculating  the  viscous/inviscid  interaction  is  to  calculate 
the  inviscid  flow  over  the  basic  body.  The  resulting  distribu¬ 
tion  of  the  velocity  at  the  boundary  is  then  prescribed  in  the 
second  iteration  step  as  the  boundary  layer  edge  velocity  ue  for 
the  boundary  layer  calculation  in  the  direct  mode  up  to  the 
separation  point.  The  boundary  layer  calculation  is  then 
switched  to  the  inverse  mode  and  the  solution  is  carried  on  into 
the  separated  flow  using  a  prescribed  distribution  of  the 
boundary  layer  displacement  thickness  5*.  The  result  of  the 
inverse  calculation  is  a  solution  for  the  boundary  layer  edge 
velocity  (the  "viscous  velocity”)  which,  in  general,  will  not 
agree  with  the  "inviscid  velocity"  produced  by  the  inviscid  flow 
theory.  The  effective  displacement  surface  between  separation 
and  a  point  downstream  of  reattachment  or  in  the  plume  entrain¬ 
ment  region  is  assumed  to  be  conical  (see  fig.  1) .  An  iterative 
procedure  is  used  to  find  the  particular  cone  angle  8g  and  the 
separation  point  location  xs  for  which  the  "viscous  velocity" 
and  the  "inviscid  velocity"  will  agree  within  an  acceptable 
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tolerance.  In  the  remainder  of  this  report,  this  procedure  will 
be  referred  to  as  the  "xs-0s  iteration". 

In  the  improved  interaction  procedure  the  separation  point 
location  xs  is  assumed  known.  Also,  in  order  to  obtain  a  first 
approximation  for  the  shape  of  the  displacement  surface  an 
angle  0s  is  assumed  known.  These  quantities  can  be  obtained 
from  the  method  of  reference  5  or  several  other  available 
methods  such  as  those  evaluated  by  Abeyounis  (ref.  9) .  The 
angle  Gs  is  usually  estimated  to  be  10°  as  a  first  approxima¬ 
tion.  Thus,  the  new  interaction  procedure  begins  in  the  same 
manner  as  in  reference  5.  The  first  step  is  to  calculate  the 
inviscid  flow  over  the  basic  body.  The  second  step  is  to  cal¬ 
culate  a  boundary  layer  displacement  thickness  assuming  a 
conical  displacement  surface  in  the  separated  region.  The  next 
step  in  the  new  interaction  procedure  is  the  computation  of  the 
new  displacement  thickness  distribution  in  which  the  assumption 
of  a  conical  surface  is  dropped  and  the  new  displacement  thick¬ 
ness  is  deduced  from  the  mismatch  of  the  viscous  and  inviscid 
velocities  by 


r*  -*  v  , 

new  old  u 

I 

The  displacement  surface  is  prescribed  in  this  manner  for  all 
subsequent  iterations.  Also,  it  is  no  longer  necessary  to 
operate  the  boundary  layer  method  in  a  direct  mode  up  to  the 
separation  point  but  the  inverse  method  is  used  with  the  dis¬ 
placement  surface  prescribed  over  the  entire  afterbody  flow 
region.  Carter  (ref.  10)  has  shown  that  this  update  procedure 
for  the  displacement  thickness  follows  logically  from  an 
approximate  analysis  of  the  von  Karman  momentum  integral  equa¬ 
tion. 


It  is  noted  that  the  extension  of  the  calculations  across 
both  the  separation  point  and  the  reattachment  point  in  the 
present  method  is  accomplished  ir.  the  same  manner  as  for  the 
method  of  reference  5.  The  approach  employs  an  engineering 
approximation  based  on  the  assumptions  that  upstream  influences 
are  transmitted  predominantly  through  the  external  inviscid  flow 
and  that  violations  of  the  boundary  layer  assumptions  such  as 
significant  normal  pressure  gradients  associated  with  separation 
can  be  neglected  on  the  grounds  that  they  are  only  important  in 
a  small  neighborhood  of  the  separation  point  and  have  negligible 
effect  on  the  rest  of  the  flow.  In  keeping  with  these  assump¬ 
tions,  the  location  of  the  separation  point  is  obtained 
independently  from  the  boundary  layer  calculation  and  the  skin 
friction  is  specified  to  be  zero  at  that  point  keeping  the 
boundary  layer  edge  velocity  and  the  displacement  thickness  6 
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continuous.  By  specifying  the  displacement  thickness  over  the 
entire  calculative  region  a  smoothly  converging  iteration  is 
produced.  The  update  procedure  on  the  displacement  thickness, 
including  a  relaxation  parameter,  w,  can  be  written  as 


f 

6*  = 

1  +  to 

ue 

v  -  i 

n 

i 

n-l 

u 

1 

e  t 

l  1  J 

where,  if  uj  =  1,  equation  (1)  is  obtained.  It  was  found  that 
overrelaxation  (oj  >  1)  could  be  used  to  accelerate  the  con¬ 
vergence  of  the  interaction  process.  However,  with  over relaxa¬ 
tion  an  oscillatory  pattern  was  found  to  be  produced  in  the 
convergence  so  that  for  practical  calculations,  underrelaxation 
is  usually  used  (u  <  1) .  In  the  remainder  of  this  report,  this 
procedure  will  be  referred  to  as  the  "6*-update  procedure". 

It  is  noted  that  the  new  procedure  does  not  provide  any 
new  information  for  locating  the  separation  point,  xg.  Abeyounis 
(ref.  8)  examined  eight  empirical  and  semi-empirical  separation 
criteria  and  prediction  methods  (not  including  that  of  reference 
5,  the  "xs-8s  iteration”)  using  experimental  data  from  boattailed 
afterbodies  with  solid  exhaust  plume  simulators.  Abeyounis 
found  that  none  of  the  methods  accurately  predicted  the  separa¬ 
tion  locations  for  all  the  conditions  included  in  the  data. 

The  method  of  reference  5  has  been  found  to  compare  favorably 
with  the  best  of  the  criteria  examined  by  Abeyounis.  The  xs-Gs 
iteration  procedure  remains  an  option  in  the  revised  computer 
code . 


3.  SUPERSONIC  INTERACTIONS 


The  interaction  between  a  boundary  layer  and  a  shock  wave 
of  sufficient  strength  to  cause  separation  of  the  boundary  layer 
is  very  complex  and  presents  formidable  difficulties  for  develop¬ 
ing  approximate  models  for  engineering  calculations.  The  most 
detailed  experimental  studies  of  shock-wave-boundary-layer  inter¬ 
actions  in  the  Mach  number  range  1.2-1. 5  have  concerned  normal 
shock  waves.  While  the  shock  waves  encountered  on  axisymmetric 
nozzle  boattails  are  not  necessarily  normal  to  the  surface,  a 
similar  degree  of  complexity  of  the  flow  to  that  associated  with 
a  normal  shock  is  to  be  expected.  Therefore,  some  understanding 
of  the  nature  of  the  flow  can  be  gained  by  examination  of  a 
normal-shock-boundary  layer  interaction.  An  important  work  for 
this  purpose  is  that  of  Seddon  (ref.  11)  which  discusses  the 
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fundamental  nature  of  normal  shock  and  turbulent  boundary 
layer  interaction  from  detailed  examination  of  velocity  profiles. 

An  interesting  feature  of  Seddon's  work  is  evidence  obtained 
that  the  wake  hypothesis  of  Coles  holds  true  over  the  greater 
part  of  the  flow  except  in  the  immediate  vicinity  of  the  shock 
wave  or  other  regions  of  strong  pressure  gradient.  This  is  an 
important  observation  for  the  present  work.  The  boundary  layer 
method  used  previously  for  subsonic  flows  is  based  on  the  Coles 
law  of  the  wall  and  law  of  the  wake  so  that  the  work  of  Seddon 
suggests  the  possibility  of  extending  the  method  to  flows  with 
shock  waves. 

In  developing  the  calculative  technique  to  be  described 
here,  it  was  found  that  straightforward  application  of  the  cal¬ 
culative  techniques  used  in  reference  5  and  discussed  in  the 
previous  section  was  not  possible.  Calculations  using  those 
techniques  were  found  to  be  unstable.  A  new  calculative  tech¬ 
nique  was  derived  which  employs  Prandtl-Meyer  theory  and  the 
semiempirical  approach  of  Mager  (ref.  12) . 

In  order  to  produce  a  stable  calculation,  it  was  found  to 
be  necessary  to  eliminate  from  iterative  calculations  the 
boundary  layer  upstream  of  the  separation  point.  Thus,  for 
supersonic  interactions,  the  first  two  steps  of  the  calculation 
are  the  same  as  for  subsonic  interactions,  that  is,  calculation 
of  the  inviscid  flow  over  the  basic  body  followed  by  a  calcula¬ 
tion  of  the  boundary  layer.  The  separation  point  is  determined 
from  the  inviscid  pressure  distribution  by  an  approximate  method 
and  the  first  boundary  layer  calculation  is  made  assuming  a 
conical  separated  displacement  surface  as  in  the  method  of 
reference  5.  In  subsequent  iterations,  the  starting  point  for 
the  boundary  layer  calculations  is  the  separation  point.  Thus, 
the  boundary  layer  upstream  of  the  separation  point  is  calculated 
only  once  and  is  assumed  to  be  the  same  for  all  subsequent 
iterations.  Downstream  of  the  separation  point  the  theory  for 
prescribed  displacement  thickness  is  modified  by  replacing  the 
arbitarily  specified  6*  distribution  with  an  equation  employing 
the  Prandtl-Meyer  function. 


(4) 
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with 

d> 


=  u  /a  cos  <fi 
e  e  e 


(7) 


In  the  implementation  of  this  theory  the  Mach  number  M  in 
the  Prandtl-Meyer  function  is  evaluated  using  the  velocity  com¬ 
ponent  u^  as  calculated  from  the  boundary  layer  theory  for  the 
current  iteration  with  the  tangential  angle,  (f ,  evaluated  from 
equation  (3)  at  the  previous  step  of  the  current  boundary 
layer  integration. 

The  separation  point  used  in  supersonic  cases  is  not  so 
easily  predicted  as  for  subsonic  flows.  In  the  present  method, 
two  options  are  available  for  predicting  the  separation  point. 

In  one  option  the  separation  point  is  determined  from  the  in- 
viscid  solution  on  the  basic  body  at  the  first  iteration.  This 
option  depends  somewhat  on  the  arbitrary  spacing  of  the  calcula- 
tive  mesh  of  the  external  inviscid  flow  in  that  the  separation 
position  is  placed  at  the  inviscid  mesh  point  previous  to  the 
point  at  which  the  maximum  Mach  number  is  calculated  on  the 
boattail .  In  the  second  option,  the  user  of  the  method  is 
allowed  to  input  any  arbitrary  separation  point  location.  With 
the  separation  point  known,  calculation  requires  an  initial 
value  for  the  angle  of  the  displacement  surface  after  separa¬ 
tion.  This  is  determined  from  the  theory  of  Mager.  This  theory 
employs  the  following  equations: 


Pg/Pi  =  1  +  0.225  YM?/[l  +  |  (v  -  1)M?] 
(pf/pi}v  =  ps/pi  +  °-18  YM?ef/[l  +  0.275  (y  -  l)M?]j 
(Pf/Pj)  =  [2  YM^  sin2  0g  -  (Y  -  1)]/(Y  +  1) 


(8) 

(9) 

(10) 
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tan  Qf  = 


-{[* 


<Y  + 


1)  M?/  (M? 


sin 


-  1) 


1  tan 


(ID 


The  required  angle  0*-  is  determined  by  an  iterative  procedure 
as  follows.  The  value  of  is  the  maximum  Mach  number  ahead 
of  the  shock,  obtained  from  the  inviscid  solution.  Then  ps/p^ 
is  determined  from  equation  (8) ,  and  the  first  approximation  to 
Pf/p^  is  determined  from  equation  (9)  with  0f  approximated  by 


so  that 


where 


ef  j Vm-  -  1  /YM-]  [lPf/Pi>  -  l] 

Pf  Pr  ^  r  i*) 

pT  *  pf  l(1  +  +  (Ps/Pi’GJj 


G  =  -0.18  ~y M?  -  1  /[l  +  0. 275  (y  -  l)M?j 


(12) 


(13) 


(14) 


Then  0g  is  computed  from  the  oblique  shock  relation,  equation 
(10) ,  and  from  equation  (11) .  Substituting  into  equation 
(9)  and  calculating  a  new  from  equation  (10)  allows  a  second 
value  of  6f  to  be  determined  from  equation  (11) .  The  second 
value  of  0f  is  used  as  the  initial  value  of  the  turning  angle  of 
the  displacement  surface  following  Mager’s  advice  that  further 
iteration  will  not  improve  the  accuracy. 

with  the  Prandtl-Meyer  equation  replacing  the  arbitrarily 
specified  <5*  distribution,  there  are  three  equations  for  the 
three  unknowns,  ue,  <S^,  and  6*  of  the  boundary  layer  theory. 
However,  the  overall  flow  problem  is  not  closed  since  the  cal¬ 
culation  of  the  flow  in  the  limited  region  of  separation  does 
not  necessarily  satisfy  the  boundary  conditions  on  the  flow  at 
downstream  infinity.  In  order  to  provide  final  closure  of  the 
problem,  the  Prandtl-Meyer  analysis  stops  at  the  reattachment 
point  and  the  calculation  method  returns  to  the  u^  prescribed 
procedure  as  in  the  theory  of  reference  5.  In  this  procedure 
the  value  of  ue  that  was  obtained  from  the  boundary-layer- 
Prandtl-Meyer  theory  is  faired  into  the  inviscid  ue  solution 
obtained  from  the  previous  iteration.  The  result  of  this  pro¬ 
cedure  is  as  in  the  subsonic  case,  two  velocity  distributions, 
one  from  the  boundary-layer-Prandtl-Meyer  theory  and  one  from 
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the  external  inviscid-f low  theory.  In  order  to  cause  these 
velocity  distributions  to  agree  as  closely  as  possible,  another 
parameter  must  be  adjusted.  The  parameter  selected  for  this 
purpose  is  the  boundary  layer  edge  velocity  at  the  separation 
point  Ugg.  Because  of  the  approximations  inherent  in  the 
present  model,  the  velocity  at  the  separation  point  for  the 
boundary- layer  calculation  cannot  be  forced  to  be  equal  to  the 
velocity  at  that  point  in  the  inviscid  flow  for  any  arbitrary 
iteration  step.  Thus,  the  value  of  the  initial  velocity  for 
the  boundary-layer  calculation  is  iterated  until  the  mean 
squared  error  between  the  boundary-layer  velocity  and  the 
inviscid  velocity  is  minimized  as  for  the  subsonic  flows  dis¬ 
cussed  previously. 

Due  to  the  approximate  nature  of  the  model,  it  is  generally 
not  possible  to  cause  the  boundary  layer  and  inviscid  solutions 
to  agree  as  closely  as  for  subsonic  cases.  While  it  is  typi¬ 
cally  possible  in  subsonic  cases  to  achieve  agreement  to  an  rms 
error  of  less  than  one  percent,  an  error  of  four- to  five  percent 
is  more  typical  of  supersonic  cases  calculated  to  date.  How¬ 
ever,  as  will  be  shown  subsequently,  the  pressure  coefficient 
distribution  determined  from  the  boundary  layer  velocity  solu¬ 
tion  compares  very  well  with  measured  values. 


4.  EXHAUST  PLUME  ENTRAINMENT  MODEL 


The  effect  of  an  exhaust  jet  plume  on  the  flow  on  an  after¬ 
body  has  been  found  to  consist  of  two  separate  effects.  First 
is  the  effect  of  the  shape  of  the  plume.  Second  is  the  effect 
of  the  entrainment  of  external  air  into  the  plume  boundary.  The 
two  effects  oppose  each  other  since  the  entrainment  of  low  speed 
external  air  by  the  high  speed  jet  tends  to  decrease  the  effec¬ 
tive  expansion  of  the  jet  boundary. 

The  shape  effect  of  the  jet  plume  is  accounted  for  by  the 
expansion  of  the  inviscid  jet  core  flow.  The  entrainment  effect 
is  determined  by  calculating  a  jet  mixing  displacement  thickness 
and  defining  an  effective  body  shape  using  that  surface. 

The  basis  of  the  new  exhaust  plume  entrainment  model  is 
the  technique  developed  by  Peters,  et  al.  (ref.  13)  in  which  an 
integral  method  for  the  exhaust  plume  mixing  layer  is  coupled 
with  a  method  of  characteristics  inviscid  jet  calculation.  The 
modifications  to  Peters'  technique  for  the  present  calculation 
method  include  use  of  a  prescribed  external  pressure  distribu¬ 
tion,  calculation  of  integrals  in  the  integral  mixing  layer 
method  across  the  mixing  layer  only,  that  is,  omitting  the 
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external  inviscid  flow  from  the  integrals,  and  determination 
of  derivatives  of  flow  quantities  at  the  inner  boundary  of  the 
mixing  layer  by  numerical  differences  rather  than  by  solution 
of  a  set  of  simultaneous  equations.  The  boundary  of  the  jet  is 
taken  to  be  the  midpoint  of  the  mixing  layer  and  the  displace¬ 
ment  thickness  of  the  mixing  layer  is  calculated  in  a  manner 
analogous  to  that  for  a  boundary  layer  (see  fig.  2).  The 
entrainment  of  the  boundary  layer  is  thus  accounted  for  to  a 
first  approximation  by  using  the  displacement  surface  of  the 
jet  as  the  effective  boundary  upon  which  to  calculate  a  boundary 
layer.  The  boundary  layer  of  the  boattail  is  then  simply 
extended  onto  this  approximate  surface  as  though  the  jet  were  a 
solid  surface.  Entrainment  is  accounted  for  by  the  fact  that 
the  displacement  surface  of  the  mixing  layer  velocity  profile  is 
a  negative  quantity  relative  to  the  jet  boundary. 


4.1  INTEGRAL  EQUATIONS 

The  basic  equations  of  the  entrainment  in  the  coordinates 
shown  in  figure  2  are: 

Continuity 


(pur)x  +  |pr[v  -  u(rw)xjj  =  0  (15) 

Axial  Momentum 

rpuux  +  rpjv  -  u(rw)xJur  =  -rpx  +  (rp&ur)r  (16) 

The  energy  equation  is  represented  by  the  approximation 

S  =  (VTte  -  L)  =  (Tti/Tte  “  "  ue)/(ui  -  ue]  (17) 

where 

rw  =  r^  +  j  b  (see  figure  2)  (18) 

Integrating  across  the  mixing/entrainment  layer  and  using  the 
continuity  equation  to  eliminate  the  radial  velocities  yields 
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W 

j  r[p(2u  -  uw)ux  +  u(u  -  uw)px]dr 


r . 
1 


and 


+  |  IrJ  -  rj>  Cpe)x  -  r„liw6w(ur)w 
+  ri0i[Vi  -  ui(r„)x](uw  -  u±>  -  0 
r 

e 

[  r[o{2u  -  ue)ux  +  u(u  -  ue)px]dr  +  \  (r^  -  r?)  (pe). 
+  riPi  [v±  -  Vrw>x](ue  -  u±)  =  0 


(19) 


(20) 


The  results  of  the  solution  of  these  equations  are  used  to 
define  a  displacement  thickness 


ri+b 


,*  [ 

6  -j 


1  - 


pu 

u 

e  ej 


r 

r 

w 


dr 


w 


(21) 


As  a  first  approximation,  it  is  assumed  that  the  effective 
boundary  of  the  jet  is  defined  by 


r 


rw  +  6 


b  +  6* 


The  boundary  rb  is  then  used  as  an  effective  solid  boundary  on 
which  a  boundary  layer  is  computed  in  the  usual  manner . 


4.2  ENTRAINMENT  LAYER  VELOCITY  PROFILES 

The  velocity  profiles  of  the  plume  mixing  layer  are 
represented  by 
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r  -| 

r  -  r . 

1  (ui  -  ue> 

1  +  COS  IT 

L 

i 

b  jj 

(22) 


4.3  DENSITY  PROFILES 

The  mixing/entrainment  zone  density  profiles  are  determined 
in  terms  of  the  temperature  profiles  for  isentropic  flow 


P 


(1  +  me)  (S  +  1) 


where  me 


(23) 


4.4  TURBULENT  EDDY  VISCOSITY 

The  term  rw;iwgw(ur)  in  equation  (19)  represents  the 
turbulent  shear  stress  at  the  half  radius  surface.  The  factor 
Bw  is  the  eddy  viscosity  factor.  The  model  for  the  eddy  vis¬ 
cosity  used  in  this  work  is  the  same  as  that  used  by  Peters, 
et  al .  (ref.  13) . 


4.5  EQUATIONS  SOLVED 

Evaluation  of  the  integrals  of  equations  (19)  and  (20)  with 
the  profiles  of  equations  (22)  and  (23)  results  in  a  set  of  two 
ordinary  differential  equations  of  the  form 


dr. 

^  _ i. 

il  dx 


A  —  +  A 
"i2  dx  Ai3 


dPe 

dx 


=  B. 


(24) 


with  i  =  1  and  2,  where  A-j_ j  and  Bj_  are  functions  of  r^,  b,  ue, 
and  u^. 

The  equations  (24)  are  numerically  integrated  using  a 
simple  first-order  predictor-corrector  method  as  described  in 
reference  5. 

Many  terms  appearing  in  the  coefficients  of  equation  (24) 
depend  on  the  flow  conditions  at  the  inner  mixing  zone  boundary. 
In  order  to  solve  the  equations,  the  following  quantities  must 
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be  evaluated  at  r^:G^, 3u/3x,  3u/3r,  3p/3x,  9p/3r.  To  provide 
these  parameters,  the  inviscid  jet  core  flow  is  computed  using 
the  irrotational  method  of  characteristics  simultaneously  with 
the  numerical  solution  of  equations  (24)  in  a  manner  similar  to 
that  of  reference  13.  This  calculation  is  carried  out  until 
the  inner  mixing  layer  radius,  r^,  becomes  zero,  or  for  a 
specified  distance  downstream  of  the  nozzle.  Downstream  of 
that  point,  the  jet  boundary  is  assumed  to  be  cylindrical. 


5 .  RESULTS 


In  this  section,  the  theory  is  compared  with  experimental 
results  on  axisymmetric  boattail-sting  and  boattail-exhaust 
plume  configurations. 


5.1  CIRCULAR  ARC  BOATTAILS  WITH  SOLID  PLUME  SIMULATORS 

Wind  tunnel  studies  of  the  configuration  shown  in  figure  3 
were  described  in  reference  14.  The  configuration  is  a  cone- 
cylinder,  with  a  circular  arc  boattail.  Three  of  the  boattails 
are  examined  here  to  illustrate  the  capabilities  of  the  calcula- 
tive  method.  The  first  example  is  for  a  boattail  (configuration 
1  in  fig.  3)  for  which  extensive  separation  is  evident  experi¬ 
mentally.  Surface  static  pressure  distributions  were  measured 
along  the  boattail  and  along  the  sting  for  free- stream  Mach 
numbers  from  0.4  to  1.2.  Boundary-layer  transition  was  tripped 
at  x/D  =  .167.  Calculations  were  made  for  Mach  numbers  of  0.6, 
0.8,  and  0.9. 

In  figure  4  is  shown  a  typical  variation  of  the  rms  error 


~N 

r  1 

S  - 

l 

u  -  u 

e  e_ 

7k  uj 

_1 

l  V  Ij 

o 

1/2 


(25) 


between  the  viscous  and  inviscid  solutions  for  configuration  1 
in  subsonic  flow.  During  the  first  four  iterations  a  conical 
displacement  surface  was  used  to  generate  an  approximate  shape, 
and  then  the  new  procedure  [eq.  (2)]  was  used  with  w  =  .5. 

Large  oscillations  in  the  error  occur  in  the  first  four  itera¬ 
tions,  followed  by  a  smoothly  convergent  sequence.  When  the 
error  is  reduced  to  a  value  of  the  order  of  one  percent, 
oscillations  again  appear,  although  of  smaller  amplitude.  The 
error  appears  to  approach  an  asymptotic  nonzero  value.  The 
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value  of  the  minimum  error  depends  on  a  number  of  factors, 
including  the  coarseness  of  the  calculation  mesh  used  for  the 
external  inviscid  flow  and  the  smoothness  of  the  boundary.  No 
consistent  relationship  was  found  between  the  value  of  the 
minimum  error  and  the  location  of  the  separation  point. 

In  figure  5  is  shown  the  effect  of  the  new  interaction 
model  on  the  calculated  pressure  distribution  for  the  model  with 
the  circular  arc  boattail  of  configuration  1  and  cylindrical 
plume  simulator.  The  separation  point  used  for  the  calculations 
was  the  value  predicted  using  the  xs-0s  iteration  method  of 
reference  5.  The  pressure  coefficients  shown  were  calculated 
from  the  viscous  velocity,  uev.  The  new  5*  update  iteration 
procedure  using  equation  (2)  with  w  =  0.5  results  in  small 
adjustments  to  the  shape  and  in  turn  to  the  pressure  distribu¬ 
tion.  In  figure  6  is  shown  other  comparisons  between  the  new 
theory  and  data  from  reference  14  for  Mach  numbers  of  0.4,  0.6, 
0.8,  and  0.9  on  the  same  boattail  as  in  figure  5.  The  location 
of  separation  was  predicted  in  each  case  as  the  location  where 


Hi  =  **/e. 


1.2 


for  the  first  and  second  boundary  layer  iterations.  The  calcu¬ 
lated  pressure  coefficients  were  determined  from  the  average  of 
the  viscous  and  inviscid  velocity  distributions,  ue  and  uei. 

The  comparisons  for  =  0.4,  0.6,  and  0.8  are  excellent. 

The  comparison  for  M0  =  0.9  is  good  although  the  expansion  of 
the  flow  at  the  beginning  of  the  boattail  is  slightly  under¬ 
predicted.  The  difference  for  =  0.9  is  probably  due  to  the 
fact  that  a  weak  shock  occurs  which  is  almost  of  sufficient 
strength  to  cause  separation.  Attempts  to  calculate  the  flows 
using  the  theory  for  shock-induced  separation  were  unsuccessful, 
but  the  improved  subsonic  theory  predicts  the  pressure  distribu¬ 
tion  well. 

In  figure  7  is  shown  the  results  of  calculations  for  two 
cases  of  shock- induced  separation  on  the  same  body  as  for  the 
results  shown  in  figures  5  and  6.  The  separation  point  in  each 
case  was  determined  by  the  procedure  described  previously.  That 
is,  x&  was  derived  to  be  the  location  of  the  inviscid  calculative 
point  previous  to  the  point  at  which  the  maximum  Mach  number  is 
calculated  on  the  boattail  at  the  first  iteration.  The  pressure 
coefficients  were  determined  from  the  viscous  velocity  distribu¬ 
tions,  uev,  for  which  the  best  agreement  with  the  inviscid 
velocity  was  obtained.  The  results  for  a  Kach  number  of  1.15 
are  very  good  with  the  predicted  shock  location  and  plateau 
pressure  agreeing  very  closely  with  the  data.  For  a  Mach  number 
of  1.3,  the  predicted  shock  location  is  very  good  while  the 
predicted  plateau  pressure  is  in  fair  agreement  with  the  data. 


17 


AEDC-TR-80-1 9 


The  afterbody  drag  on  the  model  of  figures  5,  6,  and  7  is 
shown  in  figure  8.  The  predicted  drag  is  slightly  high  for  low 
subsonic  Mach  numbers  but  is  well  predicted  for  high  subsonic 
Mach  numbers  and  underpredicted  for  transonic  Mach  numbers.  The 
prediction  is  excellent  for  the  low  subsonic  Mach  number  of 
1.15,  but  is  slightly  low  at  the  higher  Mach  number  of  1.3. 

A  comparison  with  data  on  configuration  2  is  shown  in 
figure  9.  The  boattail  of  configuration  2  is  slightly  longer 
than  that  of  configuration  1  resulting  in  a  shorter  region  of 
separated  flow.  The  predicted  pressure  distribution  is  in 
excellent  agreement  with  the  data. 


5.2  CIRCULAR  ARC  BOATTAILS  WITH  HIGH  PRESSURE  EXHAUST  JETS 

Four  cases  with  real  exhaust  jets  are  shown  in  figure  10. 
The  configuration  1  afterbody,  was  used  with  a  sonic  nozzle 
and  jet  total  pressure  to  free-stream  static  pressure  ratios  of 
2.0,  3.0,  4.0,  and  5.0.  The  data  shown  are  for  a  free-stream 
Mach  number  of  0.8.  The  calculations  were  made  for  a  nozzle 
Mach  number  of  1.01  and  pressure  ratios  of  2.1,  3.0,  4.0,  and 
5.0,  respectively.  These  adjustments  were  necessary  for  two 
reasons.  First,  the  present  method  of  characteristics  jet 
theory  cannot  calculate  for  a  nozzle  Mach  number  of  exactly  1.0. 
Second,  the  external  static  pressure  at  the  end  of  the  boattail 
tends  to  be  slightly  higher  than  the  free-stream  pressure,  pQ, 
so  that  in  order  to  achieve  a  supersonic  flow  in  the  jet  as 
required  by  the  method  of  characteristics,  a  slightly  higher 
pressure  ratio  was  necessary  for  the  lowest  pressure  than  was 
used  in  the  experiment.  The  separation  point  locations  shown 
in  figure  9  were  determined  by  the  same  criterion  used  in  the 
previous  comparisons  shown  in  figure  6.  Good  comparisons  with 
the  experimental  data  are  calculated  for  all  four  pressure 
ratios  although  the  comparison  appears  to  improve  with  increas¬ 
ing  nozzle  pressure. 

The  predicted  afterbody  drag  for  the  body  with  an  exhaust 
jet  is  compared  with  experimental  values  in  figure  11.  The  new 
theory  appears  to  underpredict  the  drag  for  low  pressure  ratios 
and  overpredict  the  drag  for  high  pressure  ratios.  The  com¬ 
parison  could  probably  be  improved  by  using  a  different  separa¬ 
tion  location  in  the  calculations. 


5.3  SUMMARY  OF  RESULTS 

The  calculative  method  developed  previously  was  improved  by 
removing  the  restriction  that  the  separation  flow  region  have  a 
conical  shape  and  by  extending  the  Mach  number  range  to  include 
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flows  with  shock- induced  separation.  Also,  an  improved  exhaust 
jet  calculative  method  was  incorporated  based  on  a  simultaneous 
solution  of  the  method  of  characteristics  for  the  inviscid  jet 
core  flow  and  an  integral  method  for  the  jet  mixing  boundary. 

The  results  of  sample  calculations  indicate  that  for  sub¬ 
sonic  flows,  the  conical  displacement  surface  is  a  very  good 
approximation  and  relaxing  the  requirement  of  a  conical  surface 
yields  only  a  small  adjustment  to  calculated  pressure  distribu¬ 
tions  while  providing  a  substantial  increase  in  computational 
speed  and  reduction  in  cost. 

Calculations  for  supersonic  flows  indicate  the  theory 
produces  a  good  approximation  for  the  pressure  distribution  and 
consequently  for  the  afterbody  drag  for  Mach  numbers  up  to  1.3. 


6.0  COMPUTER  PROGRAM  ORGANIZATION 


In  this  section,  the  general  organization  of  the  programs 
will  be  described.  Specific  information  on  data  required  for 
input  and  data  developed  for  output  will  be  described  in  sections 
7,  8,  and  9.  Sample  Job  Control  Card  decks  for  a  typical  IBM 
370  installation  are  presented  in  section  10. 

The  overall  program  consists  of  a  mainline  program  (program 
1)  and  three  main  subprograms  each  consisting  of  several  sub¬ 
routines.  The  first  main  subprogram  is  the  inviscid-f low  program 
(program  2) .  It  is  a  modified  version  of  the  program  described 
in  reference  8 .  The  second  main  subprogram  is  the  boundary- 
layer  program  (program  3) .  It  is  based  on  the  integral  theory 
described  in  reference  5.  The  third  subprogram  is  the  inviscid 
exhaust  plume,  mixing  layer  program  (program  4) .  The  mainline 
program  controls  the  iteration  between  the  other  three  programs. 
Either  of  the  first  two  subprograms  may  be  used  separately, 
without  iterating  by  appropriate  choice  of  the  input  parameters. 

Both  the  viscous-flow  programs  and  the  inviscid-f low  pro¬ 
grams  require  some  punched  card  input  and  may  require  some  input 
data  from  disc  or  tape  data  files.  The  data  files  must  be 
identified  by  specific  Logical  Unit  numbers.  The  general  rela¬ 
tionship  of  the  programs  and  the  various  data  files  are  shown  in 
figure  12.  The  specific  Logical  Unit  numbers  required  for  input 
and  output  are  listed  in  Table  I.  Two  Logical  Unit  numbers  are 
associated  with  each  data  file  shown  in  figure  12.  One  unit  is 
used  for  input,  the  other  for  output. 

Program  2  requires  initially  data  from  cards  describing 
the  free-stream  conditions,  the  computational  mesh  and  the  body 
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shape.  Alternately,  the  program  can  accept  the  input  body  shape 
from  data  file  1.  It  can  also  accept  an  initial  solution  for 
the  perturbation  velocity  potential  from  another  data  file  (data 
file  2) .  Program  2  produces  printed  output  lists  of  the  appro¬ 
priate  flow  field  quantities,  quantities  describing  the 
configurations  and  the  computational  mesh  and  several  data  files. 
Data  file  2  is  rewritten  using  the  new  solution  for  the  poten¬ 
tial.  A  third  file  (data  file  3)  is  written  containing  the 
distribution  of  the  axial  velocity  component  at  the  inviscid 
boundary  for  use  by  the  viscous  flow  program. 

Program  3  requires  initially  the  free-stream  conditions  and 
gas  constants  as  well  as  parameters  describing  the  shape  of  the 
surface  over  which  the  boundary  layer  is  flowing.  On  the  first 
iteration  of  a  viscous-inviscid  interaction,  the  surface  shape 
is  the  same  as  that  for  program  2.  On  subsequent  iterations, 
the  body  shape  for  the  boundary-layer  program  remains  the  same, 
except  when  it  is  modified  by  the  calculation  of  a  new  plume 
shape,  while  that  for  the  inviscid  flow  is  modified  by  the 
addition  of  the  boundary  layer.  Program  3  can  also  accept  data 
from  data  files  as  optional  input.  The  boundary- layer-edge 
velocity  distribution,  ue,  can  be  input  from  data  file  3,  pro¬ 
duced  by  program  2.  The  distribution  of  the  body  shape  augmented 
by  the  displacement  thickness  can  be  input  from  data  file  4. 

That  file  differs  from  data  file  1  because  it  contains  the  raw 
data  for  6*  +  rw  versus  x  as  calculated  by  the  boundary-layer 
program  while  data  file  1  contains  the  shape  interpolated  to  the 
x  stations  of  the  original  input  shape.  If  program  3  is  being 
used  separately  from  program  2  (i.e.,  without  iterating),  two 
additional  options  are  available  for  card  input.  Either  the 
boundary-layer-edge  velocity,  ue,  can  be  input  as  mentioned 
previously,  or  the  displacement  thickness,  <5*,  may  be  input. 

These  options  are  described  more  fully  in  sections  8.1  and  8.2. 

Program  3  produces  as  output  lists  of  the  boundary  layer 
and  flow  quantities  as  they  are  calculated  along  the  body.  In 
addition,  program  3  produces  an  updated  version  of  data  file  4 
and  the  augmented  body  shape  (data  file  1)  required  by  program 
2.  Data  file  4  also  contains  a  list  of  the  velocity  ratio 
ue/ueQ  corresponding  to  the  boundary  layer. 

Program  3  can  be  used  in  a  two-dimensional  boundary  layer 
mode  if  desired.  However,  calculation  of  a  viscous-inviscid 
interaction  or  exhaust  plume  entrainment  can  only  be  done  for 
axi symmetric  cases  with  the  present  inviscid  program. 

All  card  input  pertaining  to  programs  2  and  3  is  input 
through  an  input  subroutine  called  by  program  1.  Program  1  also 
produces  a  file  of  the  quantities  needed  to  restart  and  continue 
the  calculation  if  the  calculation  should  terminate  before  all 
iterations  are  completed.  These  data  are  stored  on  data  file  5. 
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Detailed  instructions  regarding  restarting  are  presented  in 
sections  10.2.3  and  11.2. 


7.  INPUT  TO  THE  PROGRAMS 


The  data  required  by  the  programs  generally  fall  into  three 
categories:  (1)  geometrical  data;  (2)  flow  field  data;  and  (3) 

control  parameters.  The  control  parameters  are  indices  for 
specifying  options  and  iteration  counters.  It  will  be  noted  by 
comparison  with  reference  8  that  a  number  of  input  quantities 
required  for  the  inviscid-f low  program  have  been  eliminated  in 
the  present  version.  This  has  been  done  by  incorporating  the 
calculations  required  to  obtain  some  of  the  quantities  into  the 
present  code  or  by  simply  defining  fixed  values  which  have  been 
found  to  be  successful.  Specifically,  a  value  of  1.4  is  used 
for  the  initial  value  of  the  subsonic  relaxation  factor,  a  value 
of  0.1  is  used  for  the  initial  value  of  the  supersonic  relaxa¬ 
tion  factor  and  a  value  of  1.3  is  used  for  the  exponent  in  the 
normal  coordinate  stretching  function.  Also,  it  is  assumed  that 
the  computational  grid  has  equal  step  sizes  in  both  coordinate 
directions  at  the  nose  of  the  body.  Finally,  the  number  of 
relaxation  steps  allowed  in  the  inviscid  program  is  fixed  at  20 
for  the  first  four  interaction  iterations  and  is  changed  to  40 
and  80  after  four  and  eight  iterations,  respectively. 

The  general  requirement  of  the  input  data  is  that  the 
tabular  lists  of  the  various  distributions  required  represent 
smooth  curves.  This  is  especially  true  of  the  list  of  body 
shape  coordinates.  The  inviscid  program  uses  cubic  splines  to 
fit  the  input  coordinates,  so  those  coordinates  must  accurately 
represent  a  smooth  curve  with  continuous  second  derivatives. 


7.1  TABULAR  FORM 

The  input  data  required  for  calculating  transonic  viscous, 
inviscid-f low  interactions  consist  of  several  punched  cards  con¬ 
taining  parameters  describing  the  free-stream  flow  conditions, 
the  computational  mesh  for  the  inviscid  calculation,  initial 
values  for  the  viscous  flow  calculation,  and  certain  options 
that  are  available  in  the  program.  A  dictionary  of  the  input 
data  is  presented  in  the  next  section.  Table  II  shows  the  input 
variables  as  they  are  to  be  punched  on  the  data  cards.  More 
detailed  explanation  of  the  requirements  for  the  inviscid-f low 
program  (program  2)  are  presented  in  reference  8  and  are  not 
repeated  herein. 
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7.2  DICTIONARY  OF  INPUT  VARIABLES 

The  variables  required  for  input  on  punched  cards  are 
defined  in  this  section  in  the  order  in  which  they  are  required. 
Additional  details  on  the  format  of  the  punched  data  are  given 
in  Table  II.  The  first  three  cards  of  any  input  data  deck  con¬ 
tain  a  description  of  the  case  being  calculated.  Any  or  all  of 
these  three  cards  may  be  blank,  but  all  three  are  required.  The 
remaining  variables  in  Table  II  are  as  follows: 

NRSTRT  Integer  indicating  whether  calculation  is  being  re¬ 
started  to  continue  a  previous  calculation.  Only 
valid  for  interaction  calculations  (LPROG  =  0,  and 
LITER  =  TRUE.  See  below). 

=  0  Start  from  zero.  Input  all  quantities  on  cards 
or  data  files  as  required. 

>  0  Restart.  Input  data  file  5  (Logical  Unit  11) 
containing  data  from  previous  iteration  plus 
all  other  input  data  files  plus  other  card  data 
as  indicated  in  Table  II. 

N3  Integer  indicating  whether  restart  is  to  begin  with 

new  values  of  x&  and  0g  (see  section  11.2  and  Table  III)  . 

ILIM  Integer  number  of  seconds  corresponding  to  estimated 

length  of  calculation.  When  the  time  from  the  begin¬ 
ning  of  the  calculation  is  within  10  seconds  of  this 
limit,  the  final  restart  files  will  be  written.  This 
assures  that  the  job  will  not  terminate  while  writing 
such  files. 

LPROG  Integer  indicating  level  of  calculation. 

=  -1  Inviscid  flow  only. 

=  0  Viscous-Inviscid  interaction. 

=  1  Boundary  layer  only. 

Nl  Integer  iteration  counter  for  inner  viscous-inviscid 

iteration 

N2  Integer  iteration  counter  for  xs-6s  iteration  presently 

limited  to  a  maximum  of  20. 


IBL 


Integer  indicating  how  interaction  calculations  are  to 
begin . 
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=  3  Start  with  inviscid  flow. 

=  0  Start  with  boundary  layer. 

IUNIT  Integer  indicating  which  value  of  the  gas  constant, 

RGAS,  the  energy  conversion  factor,  Xj,  the  accelera¬ 
tion  of  gravity,  GC,  and  the  constants  in  Sutherlands 
temperature-viscosity  relation  are  to  be  used.  The 
choice  depends  on  whether  air  is  the  gas  being  calcu¬ 
lated  and  the  units  of  the  input  quantities. 

=  1  Input  units  must  be  pounds,  feet,  seconds,  and 
°R. 

=  2  Input  units  must  be  pounds,  inches,  seconds, 
and  “R. 

=  3  Input  units  must  be  newtons,  meters,  seconds, 
and  °K. 

For  another  gas,  or  other  units,  put  in  anything  for 
IUNIT  and  put  in  nonzero  values  of  Vise,  RGAS,  XJ,  GC, 
and  SCON.  Any  units  are  allowed.  The  basic  rule  is 
that  all  input  quantities  must  be  consistent  with 
regard  to  units. 

MIT  Number  of  relaxation  cycles  allowed  for  the  inviscid- 

flow  program  (program  2) .  If  a  value  of  zero  is  input, 
a  value  of  20  is  used.  Otherwise,  the  input  value  is 
used . 

GAM  Ratio  of  specific  heats  of  the  external  inviscid  flow. 

AMINF  Free-stream  Mach  number. 

IXY  Integer  number  of  values  of  coordinate  pairs,  X0,Y0, 

to  be  input  for  inviscid  body  shape.  If  IXY  =  0,  the 
required  shape  must  be  input  from  data  file  1  (Logical 
Unit  14).  Maximum  value  is  200. 

XO , YO  Axial  and  radial  coordinates  of  body  shape  for 

inviscid-f low  calculation,  2  per  card. 

The  next  series  of  variables,  items  7  and  8  in  Table  II,  are  for 
the  inviscid-f low  program  (program  2) .  Detailed  information  on 
how  these  data  are  to  be  obtained  is  contained  in  reference  1. 

IMAX  Number  of  grid  lines  in  the  tangential  direction;  1=1 

is  the  forward  stagnation  line,  I  =  IMAX  is  the  rear 
stagnation  line  for  closed  bodies  and  downstream 
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infinity  for  open  bodies.  For  each  grid  refinement 
IMAX  is  increased  such  that  IMAXnew  =  2 (IMAXqld)  "  1- 
The  present  limit  on  IMAX  is  81.  Instructions  for 
changing  this  limit  appear  as  comments  in  the  program 
listing  (subroutine  ONEO) .  The  grid  refinement  option 
is  generally  not  used  for  viscous-inviscid  interaction 
calculations . 

JMAX  Number  of  grid  lines  in  the  normal  direction;  J  =  1 

Corresponds  to  an  infinite  distance  from  the  body  and 
J  =  JMAX  is  on  the  body.  The  same  formula  and  limit 
that  apply  to  IMAX  also  apply  to  JMAX. 

MHALF  Number  of  grid  refinements  to  be  done.  For  inter¬ 

action  calculations  a  value  of  zero  should  be  used 
with  a  grid  fine  enough  for  adequate  resolution. 

KLOSE  Body  type . 

=  0  Open  body  (i.e. ,  one  with  a  sting  or  wake) . 

=  1  Closed  body. 

LREADP  Integer  indicating  whether  initial  estimate  of  poten¬ 
tial  distribution  is  to  be  input  from  data  file  2 
(Logical  Unit  13) . 

=  0  No. 

=  1  Yes. 

If  any  one,  or  all  of  the  next  four  input  quantities  are  input 

as  zero,  the  pre-programmed  values  are  used. 

DNDZO  Step  size  of  the  normal  coordinate  at  the  body.  The 
pre-programmed  value  is  7  percent  of  the  maximum  body 
diameter. 

XIXM  Value  of  the  computational  coordinate,  X,  at  the 

matching  point  of  the  two  stretching  functions  used  in 
the  finite-difference  scheme  (see  ref.  3) ,  for  open 
bodies  only.  Since  X  varies  from  zero  to  one,  XIXM 
is  the  fraction  of  the  total  number  of  grid  points 
which  will  be  in  the  first  stretching  region  (ahead  of 
Xjjj)  .  The  pre-programmed  value  is  0.75. 

XM  Axial  location,  xm  (in  physical  coordinates) ,  of  the 

junction  (or  matching  point)  between  the  two  tangential 
stretching  functions,  for  open  bodies  only,  see  refer¬ 
ence  8.  Must  be  less  than  XO(IXY).  This  parameter  is 
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DSDXIM 

XBT 

DMAX 

XZNEW 


PLUMIN 

GAMAP 

PTPPFS 

AMP 

THETAP 

TTP 

GMP 

XJ 


used  to  concentrate  computational  mesh  points  in  a 
certain  region.  The  usual  approach  for  interaction 
calculations  (the  preprogrammed  value)  is  to  let  xra 
be  equal  to  the  length  of  the  body  to  the  beginning 
of  the  sting,  or  plume,  XBT. 

Step  size  of  the  tangential  coordinate  at  the  junction 
between  the  two  tangential  stretching  functions.  The 
preprogrammed  value  is  8  percent  of  the  length  of  the 
afterbody  (XBT-XZNEW) . 

Length  of  the  body. 

Maximum  body  diameter. 

Length  of  forebody.  Also,  this  is  the  axial  location 
at  which  boundary-layer  calculations  will  begin  after 
four  iterations  have  been  calculated.  The  usual  pro¬ 
cedure  is  to  start  the  boundary-layer  calculations 
close  to  the  nose  of  a  long  body  at  XZ  (see  item  17  in 
Table  II)  and  then  after  four  iterations  move  the 
starting  point  to  XZNEW  (XZNEW  >  XZ) .  In  subsequent 
iterations,  the  boundary  layer  does  not  change  for 
X  <  XZNEW,  For  long  slender  bodies  with  boattails, 
XZNEW  should  be  the  beginning  of  the  boattail. 

Logical  variable  indicating  whether  an  exhaust  plume  is 
being  calculated. 

Ratio  of  specific  heats  of  exhaust  gas. 

Ratio  of  nozzle  total  pressure  to  free-stream  static 
pressure . 

Nozzle  exit  Mach  number. 

Nozzle  exit  divergence  angle. 

Nozzle  total  temperature. 

Ratio  of  gas  constant  of  air  to  that  of  the  exhaust 
gas. 

Energy  conversion  factor.  If  one  of  the  preprogrammed 
values  is  acceptable,  put  in  a  value  of  0.0.  The 
value  used  will  then  be  determined  by  the  value  of 
IUNIT  on  item  number  3  as  follows: 
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XJ 

778.0 
9336.0 
1.0 

If  one  of  the  preprogrammed  values 
is  acceptable,  put  in  a  value  of  0.0.  The  value  used 
will  then  be  determined  by  the  value  of  IUNIT  on  item 
number  3  as  follows. 

IUNIT  GC 

1  32.2 

2  386.4 

3  9.81 

The  remaining  input  variables  are  related  to  the  boundary— layer 

program  alone  or  to  the  viscous-inviscid  interaction  method. 

LSEP  Logical  variable  indicating  whether  the  location  of 

separation  (XSEP)  is  known  (TRUE)  or  not  (FALSE) . 

LITER  Logical  variable  indicating  whether  the  conical  sepa¬ 
rated  interaction  procedure  is  to  be  used.  If  LITER  = 
TRUE,  the  calculation  procedure  is  the  same  as  that 
described  in  reference  5  for  bodies  with  solid  plume 
simulators  until  a  solution  is  reached.  When  the  xs- 
6S  iteration  terminates,  the  calculation  then  begins 
the  5  -update  procedure  described  in  section  2  of  this 
report.  If  LITER  =  FALSE,  the  calculation  follows 
the  procedure  described  in  section  2  or  3  of  this 
report,  depending  on  the  value  of  the  next  input 
variable,  LNSHK. 

LNSHK  Logical  variable  indicating  whether  the  shock-induced 
separation  procedure  is  to  be  used.  If  LNSHK  =  TRUE, 
the  calculation  procedure  for  subsonic  flow  (no  shock) 
is  used  even  if  supersonic  flow  with  a  shock  wave 
occurs.  This  parameter  is  used  to  override  the  shock- 
separated  flow  procedure  for  transonic  flows  where 
weak  shock  waves  occur  and  the  supersonic  procedure 
does  not  yield  a  good  result.  If  LNSHK  =  FALSE,  the 
supersonic  procedure  will  be  followed  automatically 
if  a  shock  wave  with  Mach  number  greater  than  1.25  is 
found  on  the  boattail  on  the  first  iteration. 

XSEP  Axial  location  of  separation.  Put  in  only  if  LSEP  = 

TRUE . 


IUNIT 

1 

2 

3 

GC  Gravity  constant. 
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DTHET 


IOPT 


K 


LVAR1 


LSHAPE 


LIC 


LDSTAR 


Angle  of  5*  surface  with  the  boattail  tangent  at  xs. 
Put  in  only  if  LSEP  =  TRUE.  A  value  of  0.0  will  cause 
the  program  to  use  value  of  10.0  (degrees). 

Integer  indicating  the  mode  of  the  calculation. 

=  1  ue  is  to  be  input. 

=2  6*  is  to  be  input. 

Put  in  a  value  of  1  for  starting  an  interaction  calcu¬ 
lation. 

Integer  indicating  whether  flow  is  axisymmetric . 

=  0  Two  dimensional. 

=  1  Axisymmetric. 

Integer  indicator  for  method  of  input  of  ue  when 
IOPT  =1. 

=  0  Input  u0  (dimensional)  on  cards. 

=  1  Input  ue/ueo  on  cards. 

=  2  Input  ue/ueo  from  data  file  3  on  Logical  Unit 

12* 

Integer  indicating  option  for  calculating  all  initial 
conditions  (I.C.)  except  ug  (see  section  8.4). 

=  0  Input  initial  values  per  LIC. 

=  1  Calculate  I.C.  for  flat  plate. 

=  2  Calculate  I.C.  for  cylinder. 

=  3  Calculate  I.C.  for  cone. 

Integer  indicating  initial  condition  options  for 
IOPT  =  1  and  LSHAPE  =  0. 

=  1  Put  in  CFC1  and  DELTAl . 

=  2  Put  in  CFC1  and  DELST1. 

Integer  indicating  whether  a  file  of  6*  +  rw  is  to  be 
input. 
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=0  No  input. 

=  1  File  of  6*  +  rw  versus  x  is  required  on  Logical 

Unit  15  [data  file  4). 

LSHPBL  Integer  indicating  whether  body  shape  is  to  be  input 
for  boundary-layer  calculations. 

=  0  XRL  and  RL  are  assumed  to  be  the  same  as  XO  and 
YO.  This  is  usually  the  case  when  starting  an 
interaction  calculation. 

=  1  XRL  and  RL  will  be  required. 

NOTE:  The  next  three  variables,  NVAR,  XVAR,  VAR,  are  only 
required  on  cards  if  LPROG  =  1  and  LVARl  =  0  or  1. 


NVAR 


Integer  indicating  the  number  of  values  to  be  input  for 
the  prescribed  variable  (ue  or 
100. 


6*)  . 


Maximum  value  is 


XVAR, VAR  Axial  location  and  value  of  prescribed  variable  as 
follows : 


EL 


PT 

TT 

TWONTT 

vise 


IOPT  =  1  and  LVARl  =  0,  VAR  =  ue 
IOPT  =  1,  and  LVARl  =  1,  VAR  =  ue/ueQ 
IOPT  =  2,  VAR  =  <5* 

Reference  length.  Needed  if  input  data  lengths  are 
nondimensionalized.  If  lengths  are  dimensional,  put 
in  EL  =  1.0. 

Total  pressure,  pt(lb/ft2),  (lb/in2),  or  (newton/m2). 
Total  temperature,  (°R)  or  (°K). 

Ratio  of  body  surface  temperature  to  total  temperature, 
Tw/Tt- 


Constant  A  in  Sutherlands  formula  for  viscosity. 


V 


t3/2 
T  +  T 

s 


If  one  of  the  programmed  values  is  acceptable,  put  in 
a  value  of  0.0.  The  value  used  will  then  be  deter¬ 
mined  by  the  value  of  IUNIT  on  item  number  3  as 
follows : 
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RGAS 


SCON 


DFACT 


XZ 


IUNIT  VISC 

1  2.27  ( 1 0 — 8 )  lb  sec/ft2  (°R)1/2 

2  1.5764 (10~10)  lb  sec/in2  (°R) 1/2 

3  1.4582  (10‘6)  Newton  sec/m2  (°K)1/2 


Gas  constant.  If  one  of  the  programmed  values  is 
acceptable,  put  in  a  value  of  0.0.  The  value  used 
will  be  determined  by  the  value  of  IUNIT  as  follows: 

IUNIT  RGAS 

1  1716.0  ft2/sec2  °R 

2  247104.0  in2/sec2  “R 

3  286.96  m2/ sec2  °K 

Constant  Ts,  in  Sutherlands  viscosity  law  (see  defini¬ 
tion  for  VISC) .  If  one  of  the  programmed  values  is 
acceptable,  put  in  a  value  of  0.0.  The  value  used 
will  be  determined  by  the  value  of  IUNIT  as  follows: 


IUNIT 


SCON 


1 

2 

3 


198.6  °R 
198.6  °R 
110.333  °K 


Relaxation  factor  for  adding  6*  to  body  in  the  xs-0s 
iteration  scheme.  Usual  value  is  1.0.  If  a  value  of 
0.0  is  put  in,  a  value  of  1.0  will  be  used. 


Axial  location  of  beginning  of  boundary-layer  calcula¬ 
tion. 


RLEN  Axial  location  of  end  of  boundary-layer  calculation. 

Usually  at  least  one  maximum  diameter  larger  than  XBT . 

XT  Axial  location  of  transition  from  laminar  to  turbulent 

boundary- layer  flow. 

DXP  Axial  interval  at  which  velocity  profiles  are  to  be 

printed.  If  a  value  of  0.0  is  put  in,  no  profiles  are 
printed.  In  any  case,  profiles  will  not  be  printed 
more  often  than  the  output  step  size  for  individual 
output  quantities  (see  section  9.1). 

HLIM  Limit  value  of  H-^  =  to  indicate  separation  as 

discussed  in  reference  5.  Input  of  a  blank  or  a  value 
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of  0.0  will  cause  a  value  of  1.2  to  be  used.  To  move 
the  value  of  xs  forward,  decrease  HLIM.  To  move  xs 
downstream,  increase  HLIM.  A  value  greater  than  4.0 
will  allow  the  calculation  to  proceed  to  a  separation 
singularity  if  one  should  occur.  The  calculation  will 
terminate  at  that  point.  See  section  5  for  a  discus¬ 
sion  of  the  use  of  HLIM  in  sample  calculations. 

CFC1  Value  of  skin-friction  coefficient  at  initial  boundary- 

layer  station  (see  section  8.4). 

DELTA1  Value  of  boundary-layer  thickness,  6  (compressible) ,  at 
initial  boundary-layer  station  (see  section  8.4). 

DELST1  Value  of  boundary-layer  displacement  thickness,  5*,  at 
initial  boundary-layer  station  (see  section  8.4). 

UE1  Value  of  boundary- layer-edge  velocity  ue  at  initial 

boundary-layer  station. 

DUEDX  Value  of  boundary-layer-edge  velocity  gradient  at 

initial  boundary-layer  station. 

NR  Integer  number  of  values  of  XRP  and  RL  to  be  input  for 

body  shape.  If  NSHPBL  =  0,  this  is  assumed  to  be  the 
same  as  IXY.  Maximum  value  is  200. 

XRP , RL  Axial  and  radial  coordinates  of  body  shape  for 

boundary-layer  calculation.  If  LSHPBL  =  0,  these  are 
assumed  to  be  the  same  as  XO  and  YO,  respectively. 

For  two-dimensional  configurations,  these  represent 
the  x  and  y  coordinates  of  a  surface  measured  from  a 
reference  plane. 


8.  PROGRAM  OPTIONS 


Several  optional  modes  of  calculation  are  available  through 
the  input  parameters.  A  description  of  the  options  and  the 
corresponding  values  of  the  pertinent  parameters  follows. 


8.1  BOUNDARY- LAYER  OPTION 

To  use  only  the  boundary-layer  program,  put  in  the  three 
card  description,  then  all  values  on  the  fourth  input  data  card, 
Item  2  in  Table  II,  should  be  zero.  Then  the  first  value  on  the 
fifth  input  data  card  should  be: 
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LPR0G=1  (see  Table  II,  Item  3) 


Of  the  remaining  variables  shown  on  item  number  3  in  Table 
II,  only  IUNIT  is  required.  The  remaining  cards  would  be  those 
corresponding  to  item  4  and  items  13  to  21  as  described  in  Table 
II.  With  this  option,  the  user  has  the  choice  of  specifying 
either  the  free-stream  velocity,  ue,  or  the  displacement  thick¬ 
ness,  5*,  through  the  variables  NVAR,  XVAR,  and  VAR  on  items  14 
and  15.  The  boundary-layer  calculation  can  be  restarted  at  any 
station  by  inputting  the  value  of  the  variables  at  that  station 
as  listed  in  the  output.  The  calculated  list  of  5*  +  rw  and 
ue/uSo  will  be  written  on  Logical  Unit  10  (data  file  4)  when  the 
calculation  terminates. 

Another  method  is  also  available  for  calculating  the  bound¬ 
ary  layer  alone.  The  boundary-layer  step  of  a  viscid-inviscid 
iteration  can  be  executed  separately.  The  appropriate  values  on 
the  fifth  card  would  be: 

LPROG  =  0 
Nl  =  1 
N2  =  21 
1BL  =  0 

IUNIT  =  1,  2,  or  3 

This  option  requires  LITER  =  TRUE  and  input  of  all  other  quanti¬ 
ties  as  though  the  iterative  sequence  were  to  be  completed.  With 
the  values  just  described,  only  the  boundary  layer  will  be  calcu¬ 
lated  and  then  the  run  will  be  terminated.  If  it  is  desired  to 
continue  the  iteration,  simply  make  N2  less  than  21.  The  free- 
stream  velocity  distribution  must  be  provided  on  Logical  Unit  12 
for  this  case.  All  other  optional  inputs  are  the  user's  choice. 


8.2  INV I SC ID-FLOW  OPTION 

To  use  only  the  inviscid-f low  program,  put  in  the  three  card 
description,  then  all  values  on  the  fourth  card  should  be  zero. 
Then  put  in  LPROG  =  -1  on  the  fifth  card  (item  number  3  in  Table 
II).  Of  the  remaining  values  on  that  card,  only  IUNIT  is 
required  and  that  quantity  is  required  only  if  PLUMIN  is  TRUE. 
After  the  first  five  cards  only  the  data  for  items  4  to  9  or  10 
as  described  in  Table  II  and  section  7.2  are  required  for  this 
option . 


8.3  VISCID-INVISCID  ITERATION  OPTION 

To  use  both  the  boundary-layer  and  the  inviscid-f low  pro¬ 
grams  iteratively,  put  in  LPROG  =  0  and  all  other  quantities  as 
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appropriate.  Such  iterations  can  be  started  with  only  the  body 
shape  and  free- stream  flow  quantities  known*  The  xs-9s  itera¬ 
tion  procedure  (LITER  =  TRUE)  may  be  restarted  to  continue  a 
prematurely  terminated  iteration.  Several  options  are  available 
to  the  user  for  restarting  an  unfinished  iteration.  See  section 
11.2  for  an  example  of  restarting.  The  simplest  option  is  to 
put  in  NRSTRT  =  1  as  the  first  value  on  the  fourth  input  data 
card  and  to  provide  the  required  input  data  files  on  Logical 
Units  11,  12,  13,  14,  and  15  (see  fig.  12  and  Table  I).  The 
only  other  data  required  for  restarting  are  the  three-card 
description  of  the  case.  The  calculation  then  proceeds  from 
where  the  previous  iteration  stopped.  Another  method  of  restart¬ 
ing,  which  is  valid  for  all  iteration  procedures,  would  be  to 
omit  the  restart  file  and  put  in  NRSTRT  =  0.  The  user  can  then 
vary  any  of  the  other  input  quantities,  using  the  data  files  or 
punched  cards  as  desired.  Note  that  the  xs-es  calculation 
terminates  when  N2  reaches  a  value  of  21.  The  value  of  N1 
increases  continuously  throughout  the  calculation  while  the 
value  of  N2  is  reset  to  1  each  time  the  xs-Ss  cycle  finds  a 
minimum  error  without  satisfying  the  convergence  criterion. 


8 . 4  BOUNDARY-LAYER  INITIAL  CONDITIONS 


Initial  values  of  boundary-layer  quantities  can  be  obtained 
in  several  ways.  The  user  can  obtain  values  of  the  skin-friction 
coefficient,  Cf,  and  either  the  boundary-layer  thickness,  5,  or 
the  displacement  thickness,  5*.  These  are  shown  on  item  18  in 
Table  II.  The  appropriate  values  LSHAPE  =  0  and  LIC  =  1  or  2 
are  punched  on  the  card  corresponding  to  item  13.  For  the  case 
when  no  other  source  of  this  information  is  available,  formulas 
have  been  included  in  the  program  based  on  the  Blasius  solution 
for  laminar  boundary  layers  and  based  on  one-seventh  power  law 
velocity  profiles  for  turbulent  flows.  These  formulas  are  only 
available  if  ue  is  being  specified  (IOPT  =  1) .  The  basic 
formulas  calculate  Cf  and  6  in  the  transform  plane  (incompress¬ 
ible,  two-dimensional) .  The  formulas  are  as  follows: 


Laminar  Flow 


(26) 


(27) 


32 


AG  DC-T  Ft  80-19 


Turbulent  Flow 


0.0592 


1-0.2 


x 


£ . 
i 


0.37  x 


-0.2 


x! 


(28) 


(29) 


These  formulas  provide  initial  values  for  boundary  layers  on 
flat  plates.  They  are  chosen  by  inputting  LSHAPE  =  1.  For 
other  geometries,  the  value  of  the  x  coordinate  is  transformed. 
Thus,  LSHAPE  =  2  chooses  the  values  for  a  circular  cylinder 
where 


x  = 


(30) 


where  xa  is  the  axial  coordinate  and  LSHAPE  =  3  chooses  the 
values  for  a  cone,  where 


x 


l  aj 


(31) 


These  formulas  have  been  found  to  be  quite  adequate  for  calcu¬ 
lating  flows  over  long  bodies.  Small  initial  errors  in  the 
calculated  boundary  layer  become  negligible  in  a  few  boundary- 
layer  thicknesses. 


8.5  OPTIONS  FOR  DETERMINING  THE  SEPARATION  POINT 

The  logical  variables  LSEP,  LNSHK,  and  LITER  and  the  input 
quantity  HLIM  are  used  to  select  a  particular  method  of  speci¬ 
fying  xs.  The  options  are  summarized  as  follows: 


33 


AEDC-TH  -80-1 9 


LSEP 

LITER 

LNSHK 

Instructions 

Calculation 

Mode 

T 

T 

- 

input  xs 

x  -e 
s  g 

F 

T 

— 

input  HLIM 

x  -e 
s  s 

T 

F 

T 

input  xg 

6*-update 

F 

F 

T 

input  HLIM 

£*-update 

T 

F 

F 

input  xs 

shock 

F 

F 

F 

xs  from  1st 

shock 

inviscid  solu¬ 
tion 

The  entire  Xg-Gg  iteration  procedure  is  a  prediction  method  for 
xs.  The  options  which  require  input  of  xs  or  HLIM  for  that  mode 
are  used  to  obtain  a  first  approximation  for  starting  the  xs-0s 
iteration.  In  addition,  those  options  can  be  used  with  the  £*- 
update  procedure  to  provide  a  fixed  value  of  xs.  Shock  separated 
cases  have  the  two  options  discussed  in  section  3. 


9 .  PROGRAM  OUTPUT 


9 . 1  STANDARD  -  OUTPUT 

The  complete  program  output  is  presented  in  the  following 
list.  Options  for  reducing  the  amount  of  printed  output  as 
discussed  in  reference  5  for  the  previous  version  of  the  program 
have  been  omitted  from  the  present  version  of  the  program. 

1.  Three-line  title  or  discription 

2.  List  of  all  values  of  integers  on  first  and  second 
data  cards 

3.  List  of  body  geometry  input 

4.  List  of  other  input  values  for  inviscid  flow 

5.  List  of  input  indices  for  boundary- layer  calculation 

6 .  List  of  body  shape  data  for  boundary  layer 

7 .  List  of  other  boundary-layer  input  quantities 
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8.  Plume  velocity  and  shape  distributions.  These  lists 
are  printed  at  the  beginning  and  after  the  4th  itera¬ 
tion. 

9.  Computed  geometric  parameters  in  normal  direction  for 
inviscid  flow 


J  -  normal  grid  index 

AN  -  normal  coordinate 

G  -  stretching  function  derivatives  (refs.  7  and  8} 

GH  -  stretching  function  derivatives  at  half  intervals 

10.  Computed  geometric  parameters  in  tangential  direction 

I  -  tangential  grid  index 

S  -  arc  length  along  reference  surface 

X  -  axial  coordinate 

Y  -  radial  coordinate 

THET  -  angle  of  reference  coordinate  surface,  6. 

For  closed  bodies,  9  is  the  same  as  the 
body  angle,  Gg.  For  open  bodies,  6  =  0g  on 
the  forebody  and  9  =  0  on  the  afterbody. 

THETB  -  body  angle,  90. 

AK  -  surface  curvature  on  closed  bodies.  For 

open  bodies  AK  is  the  surface  curvature  on 
the  forebody  and  AK  =  -(d2rw/dx2)  on  the 
afterbody. 

F  -  derivative  of  the  tangential  stretch 

function  (refs.  7  and  8) 

11.  Inviscid  relaxation  iteration  history 


IT  -  iteration  number 


DPMAX 
ID,  JD 


-  maximum  correction,  max 

ij 

-  I  and  J  location  of  DPMAX 


TT  , IT-1 

’ij  ‘  hj 


RMAX 


-  maximum  residual,  max|R.  -|,  where  R .  .  is 

ij  XJ 

the  right-hand  side  of  the  difference 

equation 
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IR,JR  -  I  and  J  location  of  RMAX 

ISUB, ISUP  -  indicates  if  maximum  residual  occurred  at 
a  subsonic  or  supersonic  point 

RAVG  -  average  value  of  the  residual 

RF1  -  relaxation  factor  for  subsonic  points 

QF3  -  relaxation  factor  for  supersonic  points 

NS  -  number  of  supersonic  points 

SEC/CY  -  time  for  iteration  cycle 

12.  Tabulated  values  of  surface  pressure  coefficient,  Cp, 
Mach  number,  and  axial  velocity  on  the  body  along  with 
a  rough  plot  of  Cp  along  the  body.  This  plot  is  dis¬ 
torted  in  the  axial  direction  because  it  is  for  equal 
spacing  in  the  computational  space.  The  asterisks 
show  the  level  of  sonic  C  .  The  inviscid  velocity 
distribution  is  used  to  calculate  C0  unless  otherwise 
indicated.  For  example,  the  label  ^ (AVERAGE  OF  INVISCID 
AND  VISCOUS  SOLUTIONS) "  indicates  that  Cp  is  derived 
from  1/2  (uev  +  uei) . 

13.  Drag  coefficients  by  trapezoidal  integration  of  the  Cp ' s 
on  the  real  body.  The  displacement  surface  is  removed 
for  calculation  of  the  drag. 

14.  Coordinates  x  and  y  of  the  sonic  line 

15.  Boundary- layer  reference  velocity,  ue  ,  unit  Reynolds 
number,  Reo/L  and  viscosity,  veQ.  This  is  only  printed 
on  the  first  step  (Kl  =  1) . 

16.  List  of  boundary-layer  quantities  with  profiles  at 
intervals  governed  by  DXP. 

X  -  axial  distance  from  the  nose,  x 

UTAU  -  friction  velocity,  u^. 

UE/UZ  -  boundary- layer-edge  velocity  ratio,  ue/u6o 
DSTPR  -  augmented  body  radius,  5*  +  rw 

HTR  -  transformed  shape  factor,  6*/S^ 
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DELTA  -  boundary- layer  thickness,  6 

DELST  -  displacement  thickness,  6* 

THETA  -  momentum  thickness,  0 

CF  -  skin  friction  coefficient,  Cf 

DELI  -  tranformed  boundary- layer  thickness,  5^ 

DUEDX  -  boundary-layer-edge  velocity  gradient, 

due/dx 

RW  -  body  radius,  rw 

17.  Quantities  showing  status  of  iteration 

XMAX  -  location  of  maximum  difference  between 

uev  and  uei 

UECHK  -  maximum  value  of  (uSv  -  ue;[)/ue  (percent) 

DRMS  -  value  of  rms  error  (percent) 

These  quantities  are  printed  each  cycle  during  the 
iteration  procedures  described  in  section  2  and  3  of 
this  report 

18.  Quantities  showing  status  of  xs-3s  iteration  when  the 
procedure  described  in  reference  5  is  being  used  to 
calculate  xs  and  6S  (LITER  =  TRUE).  The  following 
status  quantities  are  printed  each  time  the  boundary- 


layer. 

inviscid- 

■flow  iteration 

sequence  converges. 

XSEP 

-  the 

current 

value 

of 

xs 

DTHET 

-  the 

current 

value 

of 

0s 

DRMS 

-  the 

current 

value 

of 

the  rms  error  (percent) 

If  the  calculation  proceeds  to  a  normal  completion,  a 
message  will  be  printed  indicating  the  iteration  number  at  which 
the  minimum  error  occurred  and  the  value  of  other  pertinent 
parameters  for  that  iteration.  The  drag  value  and  inviscid  flow 
and  boundary-layer  data  representing  the  solution  will  be  found 
in  the  output  list  under  the  indicated  iteration  number. 


37 


AEDC-TR-SO-ig 


9.2  SPECIAL  OUTPUT  MESSAGES 

Several  special  messages  are  contained  in  the  output  to 
call  attention  to  specific  conditions  that  may  occur.  The 
messages  are  listed  in  this  section  with  instructions  about  what 
to  do  when  they  are  encountered. 

(1)  - DIVERGENCE  .  RMAX  EXCEEDS  RCHEK  , - 

This  message  is  printed  by  the  inviscid-f low  program  if  the 
relaxation  procedure  diverges.  Check  all  input  to  verify  that 
it  is  correct.  If  no  obvious  errors  appear,  the  difficulty  is 
probably  either  in  the  choice  of  parameters  for  the  computa¬ 
tional  mesh,  or  the  smoothness  of  the  data  defining  the  body 
shape . 

(2)  RFl  DECREASED  TO  _  BECAUSE  10-CYCLE  AVG  FOR  RMAX  INCREASED. 

This  message  refers  to  the  subsonic  relaxation  factor  in 
the  inviscid-f low  program.  The  initial  value  is  1.4.  The  value 
is  automatically  reduced  by  10  percent  if:  (1J  the  maximum 
residual,  averaged  over  10  cycles,  is  greater  than  that  for  the 
previous  10  cycles,  and  (2)  the  last  maximum  residual  occurred 
at  a  subsonic  point. 

{3)  QF3  INCREASED  TO  _  BECAUSE  10-CYCLE  AVG  FOR  RMAX  INCREASED. 

This  message  refers  to  the  supersonic  damping  factor  in  the 
inviscid-f low  program.  The  initial  value  is  0.1.  The  value  is 
automatically  increased  if:  (1)  the  maximum  correction,  averaged 
over  10  cycles,  is  greater  than  that  for  the  previous  10  cycles, 
and  (2)  the  last  maximum  residual  is  at  a  supersonic  point, 

(4)  INPUT  FROM  TAPE13  HAS  INCOMPATIBLE  DIMENSIONS 

This  message  is  printed  if  the  dimensions  of  the  solu¬ 

tion  read  from  Logical  Unit  13  (data  file  2  in  figure  12]  are 
not  the  same  as  the  values  of  IMAX  and  JMAX  put  on  item  number  7 
in  Table  II. 

(5)  METHOD  FOR  CALCULATING  UTAU  IN  DERIV  DOES  NOT  CONVERGE 

This  message  refers  to  the  iteration  used  to  solve  for  U_ 
when  5  is  prescribed  in  the  boundary-layer  calculation.  The  T 
only  known  cause  of  the  iteration  failing  to  converge  is  an 
error  in  the  input  data  (see  ref.  5  for  details) . 

(6)  DELTAI  HAS  BECOME  NEGATIVE 

STOP  INTEGRATION,  PRINT  PROFILE  AT  PREVIOUS  STEP 
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This  message  refers  to  the  transformed  boundary- layer 
thickness,  <5i-  The  error  condition  may  occur  due  to  the  initial 
integration  step  size  DXZ  being  too  large.  Another  possible 
cause  might  be  a  too  sudden  change  in  the  body  shape,  or  in  the 
prescribed  ue  or  6*  distribution. 

(7)  DELTA  HAS  BECOME  NEGATIVE 

STOP  INTEGRATION,  PRINT  VALUES  AT  PREVIOUS  STEP 

This  message  is  not  expected  to  occur  in  the  finished  pro¬ 
gram.  If  it  does,  check  the  input  data  carefully. 

(8)  METHOD  FOR  CALCULATING  INITIAL  VALUE  OF  DELI  DOES  NOT 
CONVERGE 

When  initial  values  of  Cf,  6,  or  6*  are  known,  the  calcula¬ 
tion  must  solve  an  integral  equation  for  the  initial  value  of 
the  tranformed  thickness,  6^.  This  is  done  by  iteration  in  a 
similar  manner  as  for  UT  described  in  message  (5) .  If  the 
iteration  does  not  converge,  it  is  usually  due  to  errors  in  the 
input  quantities. 

(9)  RESULTS  OF  ITERATION 

This  message  is  printed  at  the  end  of  each  step  of  the  xs- 
8S  cycle.  It  is  followed  by  the  current  value  of  xs,  XSEP,  the 
current  value  of  9S,  DTHET,  and  the  value  of  the  rms  error,  DRMS. 

(10)  ***FINAL  RESULTS*** 

BEST  SOLUTION  WAS  ITERATION  NO. _ 

This  message  is  printed  whenever  the  smallest  root-mean- 
squared  error  has  been  found. 

(11)  SKIN  FRICTION  HAS  BECOME  NEGATIVE  IN  AN  INCORRECT  MANNER. 
CHECK  ALL  INPUT  CAREFULLY 

This  message  will  be  printed  if  the  skin-friction  coef¬ 
ficient  changes  sign.  It  may  indicate  that  the  initial  estimate 
of  the  separation  point  location  was  too  far  downstream.  It  has 
usually  been  observed  to  occur  when  strong  shocks  are  present, 
or  when  too  few  relaxation  steps,  MIT  were  used  initially. 

(12)  ATTACHED  BOUNDARY  LAYER  SOLUTION  TERMINATED  AFTER  _ 

ITERATIONS 

This  message  will  be  printed  if  the  viscous-inviscid  inter¬ 
action  converges  and  no  separation  occurred  or  the  separation 
point  for  the  next  iteration  would  have  been  downstream  of  the 
end  of  the  body. 
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10.  PROGRAM  OPERATING  PROCEDURE 


In  this  section,  the  construction  of  card  decks  for  opera¬ 
tion  of  the  computer  programs  is  described.  First,  a  general 
description  of  the  operations  required  is  given.  Then  the 
specific  Job  Control  cards  needed  for  operation  on  an  IBM  370 
computer  are  listed.  The  same  card  decks  should  be  applicable 
at  any  370  installation  with  minor  modifications. 


10.1  GENERAL  JOB  CONTROL  SEQUENCE 

The  following  list  is  the  general  Job  Control  procedure 
that  would  be  required  to  run  the  programs  for  a  complete 
viscid-inviscid  interaction  calculation.  The  reader  is  referred 
to  figure  16  and  Table  I. 

1.  Create  partitioned  data  sets  for  restart  files  (files 
1-5  in  fig.  12) . 

2.  Define  units  2,  3,  8,  9,  and  10.  These  unit  numbers 
are  needed  for  output. 

3.  Define  units  11,  12,  13,  14,  and  15  if  NRSTRT  =  1  in 
the  input  data.  These  unit  numbers  correspond  to  the 
input  files.  They  contain  data  created  in  a  previous 
run . 

For  starting  an  initial  calculation,  the  partitioned  data 
sets  would  be  created  in  a  separate  operation.  Then,  since  no 
data  would  be  on  file,  only  units  2,  3,  8,  9,  and  10  need  to  be 
defined.  For  restarting  an  iterative  calculation,  all  data  files 
would  exist,  so  units  11  to  15  must  also  be  defined. 

To  execute  the  boundary-layer  program  alone,  unit  10  must 
be  defined  in  order  to  output  the  5*  +  rw  and  ue/ueo  list.  Unit 
12  must  be  defined  when  LVAR1  =  2,  and  unit  15  must  be  defined 
when  LDSTAR  =  1. 

To  execute  the  inviscid  program  alone,  units  2  and  8  must 
be  defined.  Unit  13  is  also  required  when  LREADP  =  1,  and  unit 
14  is  required  when  IXY  =  0. 


10.2  JOB  CONTROL  EXAMPLES 

In  this  section,  specific  examples  of  Job  Control  cards 
used  for  the  operations  discussed  previously  are  presented.  In 
the  examples,  the  computer  program  is  referred  to  as  "ITER”  with 
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the  source  code  names  "SITER"  and  the  load  module  or  binary 
version  named  "BITER" .  The  account  ID  used  in  the  examples  is 
WYL.XM.K01.  Logical  units  5  and  6  are  the  standard  input/output 
file  numbers-  It  is  not  necessary  to  specifically  define  these 
unit  numbers  in  the  JCL  deck.  An  estimated  core  memory  size  of 
256K  will  be  adequate  for  all  cases  except  when  the  restart 
procedure  is  being  used.  Then  the  estimate  should  be  increased 
to  512K. 


10.2.1  Creating  Partitioned  Data  Sets 

Partitioned  data  sets  for  use  as  input/output  disk  files 
must  be  created  before  the  normal  program  operation  can  proceed. 
The  following  procedure  is  suggested: 

Use  IBM  Utility  Program  IEFBR14 . 

Use  default  values  for  DCB  (DSORG=PO, RECFM=VS) . 

On  3330  disk,  use  SPACE  in  tracks  as  follows  {refer  to 
fig.  12  and  Table  I  for  explanation  of  file  numbers) : 

VELBOD  {File  3)  SPACE  =  (TRK , (2,1,10)) 

RESTRT  (File  5)  SPACE  =  (TRK, (10, 2, 10) ) 

PHI  (File  2)  SPACE  =  (TRK , (20 , 4 , 10) ) 

XOFILE  (File  1)  SPACE  =  (TRK , (4 , 1 , 10 ) ) 

DSFILE  (File  4)  SPACE  =  (TRK,  (G , 1 , 10} ) 

Example  of  creating  a  partitioned  data  set  called  VELBOD: 

//EXEC  PGM=IEFBR14 

//A  DD  DSN=WYL.XK. KOI. VELBOD, VOL=volume, 

//  UNIT-3330, DISP= ( ,CATLG) , 

//  SPACE=(TRK,  (2,1, 10)  ) 


10.2.2  Starting  an  Iteration  Sequence 


To  start  an  iteration  sequence,  unit  numbers,  2,  3,  8,  9, 
and  10  must  be  defined.  When  LITER  =  TRUE  in  the  input  data  set, 
the  following  sequence  of  cards  should  be  used: 
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//  EXEC  FORTGO,PROG= I TER, VOL= volume 


//  LIB= 'WYL.XM. KOI. BITER' 

//GO.FT02F0Q1  DD  DSN=WYL . XM. KOI . VELBOD (RUN1 ) , 
//  DISP— OLD 

//GO.FT03F001  DD  DSN=WYL . XM . KOI . RESTRT (RUN1 ) , 
//  DISP=OLD 


/ /GO . FT08F001  DD  DSN=WYL . XM . KOI . PHI ( RUN1 ) , 

//  DISP=OLD 

//GO . FT09F001  DD  DSN=WYL. XM . KOI . XOFILE { RUKl } , 


//  DISP=OLD 


/ /GO . FT10F001  DD  DSN=WYL . XM . KOI . DSFILE ( RUNl ) , 
//  DISP=OLD 

//GO. SYS IN  DD  * 


^  Input  data  cards 


/* 


In  the  example  presented  in  section  11.1,  LITER  =  FALSE 
and  the  restart  procedure  is  not  available.  However,  the  JCL 
cards  corresponding  to  unit  numbers  2,  3,  8,  9,  and  10  must  still 
be  defined  as  follows: 


//GO.FT02F001 

//GO.FT03F001 

//GO.FT08F001 

//GO.FT09F001 

//GO.FTlOFOOl 


DD  DUMMY 
DD  DUMMY 
DD  DUMMY 
DD  DUMMY 
DD  DUMMY 


10.2.3  Restarting  an  Iteration  Sequence 

The  specific  cards  used  to  perform  a  restart  of  the  calcula¬ 
tion  started  in  the  previous  section  (when  LITER  =  TRUE)  are: 
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//  EXEC  FORTGO, PROG=ITER, VOL=volume, 

//  LIB= 1 WYL . XH. K01BITER 1 

//GO. FT02F001  DD  DSN=WYL . XM. KOI . VELBOD (RUN2 ) , 
//  DISP=OLD 

/ /GO. FT03F001  DD  DSN=WYL . XM. KOI . RESTRT (RUN2) , 
//  DISP=OLD 

//GO.FT08F001  DD  DSN=WYL. XM . KOI . PHI { RUN2 ) , 

//  DISP=OLD 

//GO.FT09F001  DD  DSN=WYL. XM. KOI (XOFILE (RUN2 ) , 
//  DISP=OLD 

//GO.FTlOFOOl  DD  DSN=WYL . XM. KOI . DSFILE (RUN2 ) , 
//  DISP=OLD 

//GO . FTllFOOl  DD  DSN=WYL. XM. KOI . RESTRT (RUN1 ) , 
//  DISP=OLD,LABEL= {, , ,IN) 

//GO .FT12F001  DD  DSN=WYL. XM. KOI .VELBOD (RUNl) , 
//  DISP=OLD,LABEL={ , , ,IN) 

//GO . FT13F001  DD  DSN=WYL. XM. KOI . PHI (RUNl) , 

//  DISP=OLD,LABEL= (, , ,IN) 

//GO . FT14F001  DD  DSN=WYL. XM. KOI .XOFILE (RUNl) , 
//  DISP=OLD,LABEL=(, , , IN) 

//GO.FT15F001  DD  DSN=WYL. XM. KOI. DSFILE (RUNl) , 
//  DISP=OLD , LABEL= { , , , IN) 

//GO.SYSIN  DD  * 


/* 


These  cards  were  used  with  the  example  discussed  in  section  11.2. 
The  input  data  cards  required  for  restarting  are  summarized  in 
Table  III. 
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10.2.4  Executing  the  Boundary-Layer  Program  Alone 

The  specific  cards  used  to  perform  the  calculations  dis¬ 
cussed  in  section  11.3  are  listed  in  this  section.  In  the 
example  shown  here,  all  input  is  assumed  to  be  from  cards,  but 
the  output  list  of  6*  +  rw  and  ue/ueg  is  to  be  saved  on  unit  10. 
Unit  12  would  be  required  for  input  xf  LVARl  =  2,  and  unit  15 
would  be  required  if  LDSTAR  =  1.  The  cards  used  in  the  example 
in  section  11.3  are: 

//  EXEC  FORTGO,FROG= ITER, V0L= volume 
//  LIB=’SYL.XM. KOI. BITER' 

//GO.FT10FC01  DD  DSN=WYL.XM.K01 . DSFILE (RUNl) , 

//  DISP=OLD 

//GO. SYS IN  DD  * 


Input  data  cards 


/* 


10.2.5  Executing  the  Inviscid  Program  Alone 

The  cards  used  to  perform  the  calculations  discussed  in 
section  11.4  are  listed  in  this  section.  In  this  example,  the 
velocity  potential,  d> ,  is  input  from  unit  13,  and  the  new  solu¬ 
tion  for  (f»  is  output  on  unit  8.  The  calculated  velocity  on  the 
body  is  output  on  unit  2.  Input  from  unit  13  corresponds  to 
LREADP  =  1  in  the  card  input  data.  In  addition,  unit  14  would 
be  required  for  input  of  the  body  shape  if  IXY  =  0  in  the  card 
input  data.  The  specific  cards  used  in  the  example  are: 

//  EXEC  FORTGO , PROG=ITER , VOL=VOlume , 

//  LIB='WYL.XM. KOI. BITER' 

/ /GO . FT02F001  DD  DSN=WYL. XM. KOI . VELBOD (RUN2 ) , 

//  DISP=OLD 

/ /GO. FT08F001  DD  DSN=WYL. XM. KOI . PHI (RUN2 ) , 

//  DISP=OLD 

//GO. FT13F001  DD  DSN=WYL .XM. KOI . PHI (RUNl) , 

//  DISP=OLD, LABEL= (,,, IN) 

//GO.SYSIN  DD  * 
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f  Input  data  cards 

■v 


/* 


11.  NUMERICAL  EXAMPLES 


In  this  section,  several  example  calculations  are  pre¬ 
sented  to  aid  in  program  checkout.  An  example  is  presented  of  a 
complete  viscid-inviscid  interaction  using  the  6  -update  pro¬ 
cedure  described  in  section  2.  An  example  is  also  presented 
of  the  use  of  the  boundary-layer  program  alone  for  a  two-dimen¬ 
sional  geometry.  That  example  also  demonstrates  the  two  options 
for  boundary  conditions,  having  ue  specified  in  the  beginning  of 
the  calculation,  and  <5*  specified  in  the  second  part.  Input 
data  for  another  sample  case  are  also  presented  to  demonstrate 
the  use  of  the  program  to  calculate  the  inviscid  flow  alone. 


11.1  AXISYMMETRIC  INTERACTION 

A  list  of  the  punched  card  input  data  for  a  sample  calcula¬ 
tion  on  the  boattailed  body  shown  in  figure  3a  and  b  is 
presented  in  figure  13.  The  case  being  calculated  is  for  a 
free-stream  Mach  number  of  0.8.  The  body  corresponds  to  the 
ogive-cylinder  body  with  a  circular-arc  boattail  described  in 
reference  14.  The  JCL  card  deck  for  this  case  has  been  pre¬ 
sented  in  section  10.2.2.  The  running  time  for  the  complete 
calculation  is  about  3.0  minutes  on  the  IBM  370/165.  The  cal¬ 
culations  shown  were  made  by  an  IBM  3033  computer  in  1.04 
minutes.  Small  numerical  differences  will  occur  when  the  pro¬ 
gram  is  run  on  a  different  machine. 

Selected  output  for  the  sample  case  is  shown  in  figure  14. 

A  total  of  14  iterations  were  performed,  with  the  minimum  error 
of  0.74  percent  occurring  on  the  13th  iteration. 

The  complete  list  of  output  for  this  case  consisted  of  a 
total  of  approximately  5000  lines.  Output  steps  1-7  as  listed 
in  section  9.1  have  been  omitted  from  this  presentation  since 
they  simply  verify  the  input  data.  The  output  pages  shown  are 
those  corresponding  to  steps  8-18  of  the  set  described  in 
section  9.1.  Finally,  figure  14  concludes  with  the  inviscid 
solution  and  boundary- layer  solutions  corresponding  to  the  best 
result,  iteration  number  13,  and  the  final  page,  from  iteration 
14  showing  the  final  messages. 
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LI. 2  EXAMPLE  OF  RESTARTING  AN  INTERACTION  CALCULATION 

For  a  case  where  the  xs-0g  iteration  is  being  used,  the 
calculation  may  terminate  prematurely.  The  punched  card  input 
data  for  restarting  such  a  case  is  presented  in  figure  15.  The 
JCL  card  deck  for  this  calculation  was  presented  in  section 
10.2.3.  Note  that  the  iterative  calculations  cannot  be  restarted 
at  any  arbitrary  iteration  using  the  restart  file,  unit  11.  That 
file  and  the  other  output  files  contain  only  the  data  that  were 
output  just  prior  to  the  termination.  It  will  be  recalled  that 
this  restart  procedure  is  only  available  when  the  input  quantity 
LITER  =  TRUE. 


11.3  TWO-DIMENSIONAL  BOUNDARY  LAYER 

A  list  of  the  punched  card  input  data  for  a  sample  calcula¬ 
tion  on  the  two-dimensional  configuration  shown  in  figure  16  is 
presented  in  figure  17.  Note  that  two  sets  of  input  are  pre¬ 
sented  for  this  case,  giving  an  example  of  the  options  of 
prescribed  ue  and  prescribed  5*.  The  data  for  ue  and  i5*  were 
obtained  from  the  experimental  results  of  reference  15  which 
indicate  separation  occurring  in  an  adverse  pressure  gradient 
region  downstream  of  a  shock  wave.  The  output  for  the  complete 
boundary-layer  calculation  are  presented  in  figure  18.  No 
external  data  files  were  used  for  input  for  this  case.  The  JCL 
card  deck  for  this  case  was  presented  in  section  10.2.4. 


11.4  AXISYMMETRIC  INVISCID  FLOW 

The  punched  card  input  data  for  a  sample  calculation  of  the 
inviscid  flow  alone  are  presented  in  figure  19.  The  output  for 
this  case  is  shown  in  figure  20. 
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Figure  1.-  Effective  body  shape  for  separated  flow. 
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2.-  Schematic  of  exhaust  plume 
entrainment  quantities. 


50 


All  dimensions  are  referenced 
to  the  maximum  diameter 


(a)  Forebody. 


(b)  Boattail  configuration  1.  (c)  Boattail  configuration  2. 


(d)  Boattail  configuration  3. 

Figure  3.-  Body  with  circular  arc  boattail  (ref.  14). 
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Figure  5.-  Effect  of  new  interaction  model 
on  separated  flow  predictions. 
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Figure  7.-  Comparison  of  theory  with  data  for  supersonic  flows. 
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NASA  DATA  COMPARISON  -  CONFIGURATION  1 
JET  SIMULATED  WITH  HIGH  PRESSURE  AIR 
BOATTA1L  L/O  “0.8  DE/D  =0.51 


0 

0 

700 

0 

1 

1  1 

1.4 

0.8 

98 

0.0 

0.  0 

0.1 

0. 025013 

0.2 

0.  050026 

0.3 

0.075038 

0.4 

0. 100051 

0.5 

0. 125064 

0.6 

0. 150077 

0.7 

0.  175090 

0.8 

0.  200  10  2 

0.9 

0.  225115 

1.0 

0.  250128 

1.  1 

0.  275141 

1.  2 

0. 300154 

1.3 

0.  325166 

1.4 

0. 350179 

1.5 

0. 375192 

1.6 

0.400205 

1.  7 

0.425216 

1.8 

0.450230 

1.9 

0.473630 

2.0 

0.489760 

2. 1 

0.496366 

2.  2 

0. 500000 

2.  3 

0.500000 

2.5 

0.500000 

3.0 

0. 500000 

3.5 

0. 500000 

4.0 

0. 500000 

5.0 

0. 500000 

6.0 

0. 500000 

7.0 

0.  503030 

7.5 

0. 500000 

8.0 

0.500000 

B.  025 

0.499781 

6.05 

0.499125 

8.  075 

0.498030 

8.  1 

0.496496 

B.  125 

0.494521 

8.  15 

0.492104 

B.  175 

0.489241 

B.  2 

0.485931 

B.  225 

0.482171 

(a)  First  49  cards. 


Figure  13.-  Input  data  for  viscous-inviscid 
interaction  calculation. 
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S.  25 

0.477956 

8.  275 

0. 473282 

8,  3 

0.468146 

8.  325 

0. 462542 

8.  35 

0. 4564S3 

8.  375 

0. 449905 

8.4 

0. 442859 

8.  425 

0. 435319 

0.  45 

0.427277 

8.  475 

0.418722 

8.  5 

0.409646 

8.  525 

0.400037 

8.  55 

0.  389805 

8.  575 

0.  379176 

B.  6 

0.  367897 

8.  625 

0.  356032 

8.65 

0.  343565 

8,  675 

0.330479 

8.  7 

0. 316753 

8.725 

0.  30 236 B 

8.  75 

0. 287298 

8.  775 

0. 271518 

8.  8 

0. 255000 

8.  81 

0. 255000 

8.  82 

0. 255000 

8.  S3 

0. 255CC0 

8.  S4 

0. 255CCC 

3.  85 

0.  255-000 

8.  86 

0. 2550CG 

8.  87 

0. 255000 

8.88 

0.  255000 

8.89 

0. 255000 

8.9 

0. 255000 

8.92 

0. 255000 

8.94 

0. 255000 

B.  96 

0. 255000 

8.98 

0. 255000 

9.0 

0. 255000 

9.025 

0.255000 

9.05 

0.255000 

9.075 

0.  2550-00 

9.  1 

0. 255000 

9.  15 

0. 2550CC 

9.  2 

0.255000 

9.  3 

0.255000 

9.4 

0.255000 

9.  5 

0. 255000 

9.  7 

0. 255000 

9.  9 

0. 255000 

(b)  Next  49  cards. 


Figure  13.-  Continued 
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10. 1 

0. 255000 

10.  3 

0. 255000 

10.  5 

0. 255000 

io.  e 

0. 255000 

11.2 

0. 255000 

11.6 

0. 255000 

12.0 

0. 255000 

81 

31 

0  0 

0 

0.  0 

0.70 

8  .5 

T 

1.4 

2. 1 

1  .01 

F 

F 

T 

1 

1 

2  3 

1  0 

D.  5 

2134. 1 

577.8 

0. 15 

10. 0 

0. 166667 

0.  05 

8.  8 

1.0 

8.0 

0.0 

537.9 

1.0 

0.97 

O.C 

1.2 

(c)  Remaining  15  cards. 


Figure  13.-  Concluded. 
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INITIAL  VALUES  FOR  PLUME 


K 

A(  1 »K ) 

Al  2  >K  I 

A(  J,K1 

At  4pK  1 

1 

a.e ooo 

0-2550 

0.5221 

0.0 

2 

8,8000 

0.2422 

0.5221 

0.0 

I 

6.8000 

0.2295 

0.5221 

0.0 

A 

8.8000 

0. 21 57 

0.5221 

0.0 

5 

8.8000 

0.2040 

0.5221 

0.0 

6 

8.8000 

0.1912 

0.5221 

0.0 

7 

8.8000 

0.1785 

0.5221 

0.0 

8 

8  9000 

0.1457 

0.5221 

0.0 

9 

6.8000 

0. 1510 

0.5221 

0.0 
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(a)  Plume  initial  conditions. 

Figure  14.-  Selected  output  for  viscid- inviscid  interaction  calculation. 
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PLUME  CALCULATION  AT  ITERATION  NO.  1 
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(b)  Initial  plume  solution. 


Figure  14.-  Continued 
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NORMAL  COORD.  STRETCH  FOR  ALT=  1.300 
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Figure  14.-  Continued. 
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(d)  Computed  tangential  geometric  parameters 
Figure  14.-  Continued. 
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(e)  Inviscid  iteration  results  from  first  inviscid  solution. 
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(f)  First  boundary  layer  calculation  (incompl eto) . 
Figure  14.-  Continued. 


A  EDC-TR-80-1 9 


ITERATION  NO.  9  (AVERAGE  OF  I NV I  SC  ID  AND  VISCOUS  SOLUTIONS) 


PLOT  OF  CP  AT  EQUAL  XI-INCBEMENTS 


I 

XB 

TB 

FM 

qs  CP 

1 

0.0 

0.0 

0.0 

0.0  1.1704  ♦ 

ft 
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0.069 

0.017 

0.538 

0.694  0.5626 

* 

ft 
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0.137 

0.034 

0.564 

0.726  0.5091 
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* 

4 

0-207 

0.052 

0.640 

0.817  0.3499 

♦ 
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0-260 

0.070 

0.663 

0.844  0.3004 

+ 

ft 
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0.356 

0.089 

0.677 

0.860  0.2707 

4 

ft 

7 

0.435 

0.109 

0,687 

0.872  0.2495 
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ft 

8 

0.519 

0.131 

0.695 

0.881  0.2327 
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ft 

9 

0.608 

0.153 

0.701 

D.8S9  0.2177 
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ft 

10 

0.703 

0.177 

0.708 

0.896  0.2034 

4 

ft 

11 

0.804 

0.202 

0.715 

D.904  0.1866 

4 

ft 

12 

0.911 

0.229 

0.722 

0.912  0.1732 

♦ 

ft 

13 

1.026 

0.256 

0.730 

0.921  0 . 1558 

4 

ft 

14 

1.149 

0.289 

0.739 

0.931  0.1354 

4 

ft 

15 

1.280 

0.322 

0.750 

0.944  0.1103 

4 

ft 

16 

1.418 

0.357 

0.765 

0.961  0.0770 

4 

ft 

17 

1.565 

0.393 

0.768 

0.936  0.0276 

4 

ft 

16 

1.720 

0.432 

D.  838 

1.042  -0.0845 

4 

ft 

19 

1 .885 

0.471 

0.955 

1.166  -0.3389 

4  * 

20 

2.060 

0.496 

1.001 

1.213  -0.4372 

4 

21 

2.246 

0.503 

0.913 

1 . 122  -0.2489 

4  ft 

22 

2.440 

0.504 

0.652 

1.056  -0.1140 

4 

ft 

23 

2.642 

0.505 

0.832 

1.036  -0.0719 

4 

ft 

24 

2. 851 

0.505 

0.822 

1 .025  -0.0493 

4 

ft 

25 

3.067 

0.506 

0.816 

1.018  -0.0355 

4 

ft 

26 

3.200 

0.506 

0.812 

1.014  -0.0276 

4 

ft 

27 

3.516 

0.507 

0.810 

1.011  -0.0227 

4 

ft 

2S 

3.751 

0.507 

0.803 

1.009  -0.0186 

4 

ft 

29 

3.986 

0.506 

0.807 

1.008  -0.0153 

4 

ft 

30 

4.229 

0.508 

0.606 

1.006  -0.0130 

4 

ft 

31 

4.471 

0.509 

0.805 

1.006  -0.0114 

4 

ft 

32 

4.715 

0.589 

0.605 

1.005  -0.0102 

4 

ft 

33 

4.958 

0.510 

0.804 

1.005  -0.0095 

4 

ft 

34 

5.201 

0.510 

0.804 

1.005  -0.0091 

4 

ft 

35 

5,442 

0.510 

0.804 

1 .005  -0.0092 

4 

ft 

36 

5.679 

0.511 

0.604 

1.005  -0.0096 

4 

ft 

37 

5.912 

0.511 

0.805 

1 .005  -0.0104 

4 

ft 

38 

6.139 

0.512 

0.605 

1.006  -0.0116 

4 

ft 

(g)  Solution  after  4  iterations 
Figure  14.-  Continued. 
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3? 

6.359 

0.512 

0.806 

1.007 

-0.0133 

* 

ff 

40 

6.572 

0.512 

0.807 

1.008 

-0.0156 

4 

* 

41 

6.776 

0.513 

0.609 

1.010 

-0.0197 

♦ 

ff 

42 

6.®69 

0.513 

0.812 

1.013 

-0.0264 

♦ 

* 

43 

7.  152 

0.513 

0.816 

1.018 

-0.0356 

♦ 

41 

44 

7.324 

0.513 

0.821 

1.023 

-0.0455 

* 

N 

45 

7.483 

0.513 

0.827 

1.029 

-0.0593 

4 

* 

46 

7.630 

0.513 

0.83B 

1  .042 

-0.0840 

4 

* 

47 

7.763 

0.513 

0.854 

1.059 

-0.1200 

4 

ff 

46 

7.384 

0.512 

0.877 

1.084 

-0.1698 

4 

• 

49 

7.491 

0.511 

0.928 

1.138 

-0.2610 

4  * 

50 

8.086 

0.507 

0.964 

1.193 

-0.4008 

4# 

51 

8.169 

0.500 

0.997 

1  .202 

-0.4277 

4 

5? 

a. an 

0.491 

0.978 

1.176 

-0.3876 

4* 

53 

8.501 

0.480 

0.944 

1.133 

-0.3162 

4  ff 

54 

8.358 

0.468 

0.908 

1  .090 

-0.2386 

4  ff 

55 

S. 407 

0.456 

0.865 

1  .043 

-0.1455 

4 

* 

56 

8.453 

0.446 

0.821 

0.993 

-0.0464 

f 

ff 

57 

8.500 

0.433 

0.775 

0.944 

0.0543 

4 

ff 

53 

8.550 

0.413 

0.760 

0.937 

0.0695 

4 

ff 

59 

8.605 

0.412 

0.767 

0.945 

0.0722 

4 

ff 

60 

8.664 

0.401 

0.770 

0.949 

0.0667 

4 

ff 

61 

6.730 

0.358 

0.770 

0.950 

0.0654 

4 

ff 

62 

8.803 

0.374 

0.762 

0.938 

0.0648 

4 

ff 

63 

8.833 

0.358 

0.741 

0.917 

0.1298 

4 

ff 

64 

8.974 

0.340 

0.723 

0 .8*6 

0.1712 

4 

ff 

65 

9.075 

0.320 

0.690 

0.858 

0.2437 

4 

« 

66 

9.190 

0.304 

0.690 

0.872 

0.2438 

4 

ff 

67 

9.321 

0.295 

0.722 

0.911 

0.1726 

4 

ff 

68 

9.473 

0.290 

0.742 

0.934 

0.1288 

4 

ff 

69 

9.650 

0.265 

0.758 

0.952 

0.0938 

4 

« 

70 

9.859 

0.281 

0.785 

0.960 

0.0784 

4 

ff 

71 

10.110 

0.279 

0.774 

0.971 

0.0578 

4 

ff 

72 

13.417 

0.279 

0.766 

0.985 

0.0308 

4 

ff 

73 

10.800 

0.279 

0.792 

0.991 

0.0184 

4 

ff 

74 

11.293 

0.279 

0.795 

0.995 

0.0109 

4 

ff 

75 

11.950 

0.279 

0.797 

0.997 

0.0062 

4 

ff 

76 

12.870 

0.279 

0.799 

0.998 

0.0033 

4 

ff 

77 

14.250 

0.279 

0.799 

0.999 

0.0015 

4 

ff 

78 

16.550 

0.279 

0.800 

1  .000 

0.0006 

4 

« 

79 

21.150 

0.279 

0.600 

1.000 

0.0001 

4 

ff 

80 

34.948 

0.279 

0.600 

1  .000 

0.0000 

4 

« 

Al  »•««•<•••» 

0.279 

0.800 

1  .000 

0.0 

4 

ff 

(g)  Continued. 


Figure  14.-  Continued 


AEDC-TR-SO-ig 


TOTAL  BODY  DRAG  COEFFICIENTS  0.05813 


AFTERBODY  DRAG  COEFFICIENTS  0.0A53S 


STATUS  OF  ITERATION 

XHAX  =  6.4071 
UECHK  *  6.2636 
DBMS  =  2 .6226 


+- 


(g) 

Figure 


Concluded. 

14.-  Continued. 


AEOC-TR -BO-19 


INITIAL  VALUES  FOB  PLUME 


-J 


K 

At  1  ,K1 

AI2.K) 

At  3 iK 1 

A(4.K> 

1 

e.eooo 

0.2550 

0  5221 

0.0 

2 

8.6000 

0.242Z 

0.5221 

0.0 

3 

e  eooo 

0.2295 

0.5221 

0.0 

A 

6.8000 

0.2167 

0.5221 

0.0 

5 

e.eooo 

0.2040 

0.5221 

0.0 

6 

a.soao 

0.1912 

0  .5221 

0.0 

7 

8.8000 

0.1765 

0.5221 

0.0 

S 

e.eooo 

0.1557 

0.5221 

0.0 

9 

a.soao 

0.1530 

0.5221 

o.a 

10 

8.8000 

0.1402 

0.5221 

0.0 

1 1 

8.8000 

0.1275 

0.5221 

0.0 

1 2 

e.eooo 

0.1147 

0  5221 

0.0 

13 

8  6000 

0.1020 

0.5221 

0.0 

14 

e.eooo 

0.0892 

0 .5221 

0.0 

15 

8.8000 

0.0745 

0.5221 

0.0 

15 

e.eooo 

0.0557 

0.5221 

0.0 

17 

8.8000 

0.0510 

0.5221 

0.0 

13 

e.eooo 

0.0582 

0.5221 

0.0 

19 

e.eooo 

0.0255 

0  5221 

0-0 

20 

e.eooo 

0.0127 

0.5221 

0.0 

El 

e.eooo 

0.0 

0.5221 

0.0 

I 

XPEUI 

TPE(I) 

1 

6.7300 

1442.714 

2 

8.6026 

1466. 031 

3 

8.8833 

1463.626 

4 

6.9735 

1510.610 

5 

9.0750 

1579.031 

6 

9.1900 

1577.453 

7 

9,3214 

1524.848 

8 

9.4731 

1437.343 

9 

9.6500 

1462.437 

10 

9.6591 

1452.112 

1 1 

10.1100 

1436.264 

12 

10.4166 

1419.318 

13 

10.3000 

141 1 .542 

14 

11.2928 

1406.852 

15 

11.9499 

1403.899 

16 

12.8699 

1402.060 

17 

16.2499 

1400.970 

ie 

1 6 . 5498 

1400.381 

t  9 

21.1495 

1400.110 

2D 

34.9430 

1400.026 

21 

1400.018 

(h)  Plume  solution  rfter  4  iterations. 
Figure  14.-  Continued. 


AEDC-TR-80-19 


PLUIE  CALCULATION  AT  ITERATION  NO.  4 


X 

RI 

8 

RU 

UI/UE 

8.8510 

0.25436 

0.00106 

0.25480 

1.2929 

8.9020 

0 . 25368 

0.00201 

0.25452 

1 . 3005 

8.9530 

0 . 25282 

0.00304 

0.25407 

1.3114 

9.0040 

0.25180 

0,00424 

0.25352 

1 . 3426 

9.0550 

0.25068 

0.00558 

0.25289 

1.3878 

9.1060 

0.24947 

0.00705 

0.25223 

1 .4059 

9.1570 

0.24324 

0.00853 

0.25159 

1 .4050 

9,2030 

0.24703 

0.01001 

0.25093 

1 . 3943 

9.2590 

0.24536 

0. 01 142 

0.25043 

1 . 3667 

9.3099 

0.24476 

0.01273 

0.24993 

1.3405 

9.3609 

0.24371 

0.01348 

0.24944 

1 . 3228 

9.4119 

0. 24270 

0.01517 

0.24397 

1.3077 

9.4629 

0.24173 

0.01632 

0.24852 

1 .2930 

9.5139 

0.24078 

0.01742 

0.24307 

1.2837 

9.56  49 

0.2J9E6 

0.01849 

0.24762 

1 .2756 

9.6159 

0.23895 

0.01 954 

0.24719 

1 .2677 

9.6669 

0.23807 

0.02057 

0.24676 

1 .2616 

9.7179 

0.23719 

0.02158 

0.24632 

1.2589 

9.7689 

0.23632 

0.02259 

0.24539 

1 .2562 

9.B199 

0.23545 

0.02359 

0.24546 

1.2535 

9.8709 

0.23459 

0-02458 

0.24503 

1.2507 

9.9219 

0.23373 

0.02556 

0.24461 

1.2473 

9.9729 

0.23189 

0.02653 

0.24420 

1 . 2439 

10.0239 

0  23205 

0.02749 

0.24379 

1.2406 

TINE  TO  CALCULATE  PLUME  =  1.42  5ECS 


(h)  Concluded. 


Figure  14,-  Continued 


A  EDC-TR  -60-19 


ITERATION  NO.  13 


X 

UTAU 

DEIST 

UE/UZ 

DSTPR 

HTR 

DELTA 

TNETA 

CF 

DELI 

DUE  OX 

RM 

8.0Q29E+00 

3.69028*01 

t. 09228-02 

1.1422 

0.5109 

0  6889 

1 .0467E-0I 

5.6459E-03 

2.631E-03 

5.0246E-02 

4.40558+02 

0.5000 

8.04448+00 

4.001 5E+01 

1 .0420E-02 

1.1599 

0.5097 

0.8576 

1.07288-01 

5. 2556E-03 

2.889E-03 

5. D646E-02 

4.71028+02 

0.4993 

8.08598+00 

4.I076E+01 

9. 97058-03 

1.1772 

0.5073 

0.8259 

1 . I024E-01 

4.8922E-03 

2 . 939E-03 

5.11 9SE-02 

4. 60558+02 

0.4974 

6.12738+00 

4.11 37E  *  0 1 

1 .0057E-02 

1.1793 

0.5044 

0.8260 

1 .1077E-01 

4.9256E-03 

2.9358-03 

5. 1 073E-02 

7.3560E+01 

0.4943 

8.1630E+00 

4.11 87E*01 

1 .0164E-02 

1.1815 

0.5001 

0.8261 

1  . II56E-01 

4 . 9740E-03 

2.9308-03 

5.09548-02 

6.2985E+01 

0.4900 

8.20488400 

4.041 1E*01 

1  .07018-02 

1.1711 

0.4959 

0.8521 

1  . 1027E-01 

5.351 7E-03 

2,8808-03 

5.03268-02 

-2.93148+02 

0.4852 

8.24088*00 

3.9605E+01 

1.12878-02 

1.1603 

0.4908 

0.8768 

1.09538-01 

5.75188-03 

2.827E-03 

4.9694E-02 

-3.0I70E+02 

0.4795 

6.27208*00 

3.82208*01 

1 .2207E-02 

1.1410 

0.4860 

0.9159 

1 .08158-01 

6 . 3840E-03 

2.741E-03 

4.9502E-02 

-6. 2442E  +  02 

0.4738 

8.30338*00 

3.67928+01 

1 .3240E-02 

1.1210 

0.4807 

0.9540 

1 .0749E-01 

7.06938-03 

2 . 6468-03 

4.94118-02 

-6.2734E*02 

0.4674 

8.33068*00 

3. 5260E+01 

1  .4391 E-02 

1 .0999 

0 . 4756 

0.9926 

1 .0726E-01 

7.81 26E-03 

2.542E-03 

4.95738-02 

-8.9099E*02 

0.4612 

8.35808*00 

3.37028+01 

1.57118-02 

1  .0781 

0.4701 

1.0319 

1 .0759E-01 

8.6382E-03 

2.4308-03 

4.99788-02 

-9.2379E*02 

0.4544 

8.36268*00 

3, 221*8+01 

1 . 7082E-02 

1 .0576 

0.4649 

1.0688 

1 .0837E-01 

9.4702E-03 

2.320E-03 

5.0554E-D2 

-7.9508E+02 

0.4478 

8.40718*00 

3.06  7  78  *  01 

1 .65538-02 

1.0367 

0 . 4594 

1.1072 

1 .0963E-01 

1 .0393E-02 

2.203E-03 

5.1330E-02 

-7.9965E+02 

0.4407 

a.4303E*00 

2. 9205E+01 

2.03278-02 

1 .0167 

0.4540 

1.1449 

1.11 28E-01 

1 .13498-02 

2. 0068-03 

5.22468-02 

-7.9622E*02 

0.4336 

8.45348*00 

2.76998+01 

2.22538-02 

0.9965 

0.4404 

1.1848 

1 .1340E-01 

1 .241 18-02 

1.9648-03 

5.33508-02 

-7-9334E+02 

0.4261 

8.47678*00 

2.60738*01 

2.45958-02 

0.9750 

0.4427 

1  .2304 

1 . 1609E-01 

1.36558-02 

1 ,8278-03 

5.47328-02 

-8.2590E+02 

0.4181 

a. 50008*00 

2.44898+01 

2.72298-02 

0.9545 

0.4369 

1.2776 

1.19338-01 

1 .5000E-02 

1 .691E-03 

5.6286E-02 

-7.9077E+02 

0.4097 

8.52508*00 

2.32208*01 

2.9312E-02 

0.9384 

0.4299 

1 .3171 

1  -2305E-01 

1 .6285E-02 

1 .5788-03 

5.77408-02 

-6  *31 27E  +  02 

0.4000 

8.51258*00 

0,0 

2 .64328-02 

0.9384 

0.4333 

4.0000 

5 , 0946E-02 

6 . 4299E-03 

0.0 

2.1  7948-02 

-3.926  78*02 

0.4049 

8.53128*00 

-3. 14158-01 

3-05038-02 

0.9312 

0.4200 

4. 0296 

5.42Q3E-02 

6 . 0631E-O3 

-2.939E-07 

2.3057E-02 

-3.621 1E+02 

0 . 3975 

6.55008*00 

-1 . 09758*00 

3. 3537E-02 

0.9253 

0.4235 

4.1941 

5.8701 E-02 

7.2939E-03 

-3.6388-06 

2.4737E-02 

-2.9916E+02 

0.3899 

8.57732*00 

-2.3545E*00 

3.B374E-02 

0.9201 

0.4170 

4.5581 

6.5704E-02 

7.8564E-03 

-1 .6968-05 

2.7256E-02 

-1.74118*02 

0 . 3782 

8.60458*00 

-3. 51548+00 

4.5333E-02 

0.9166 

0.4111 

5.0125 

7.3882E-02 

8.4393E-03 

-3.B11E-0S 

3.01 09E-02 

-1  .2O76E  +  02 

0.3657 

8.63448+00 

-4.34148*00 

5.2579E-02 

0.9131 

0.4039 

5.4403 

8.2387E-02 

9. 1 253E-03 

-5.8628-05 

3-30428-02 

-1 . 15018+02 

0.1513 

8.66438*00 

-5. J876E+00 

6.1762E-02 

0.9136 

0.3979 

6.1165 

9. 3355E-02 

9.6869E-03 

-9.019E-0S 

3.62498-02 

4. 3254E+00 

0.3361 

8 .6971  E  *00 

-6. 00528+00 

7.21718-02 

0.9136 

0.3905 

6.7277 

1 .C489E-0I 

1 .04358-02 

-1 .15IE-04 

3.9446E-02 

-1 .36B4E+01 

0.3163 

8.7300E+00 

-6,50008*00 

6.3744E-02 

0.9136 

0.3031 

7.2833 

1 .170SE-01 

1 .13168-02 

-1 .345E-04 

4.25028-02 

-2.37918+01 

0.2994 

8 . 76638*00 

-7. 1 720E+00 

1.00108-01 

0.9149 

0.3772 

8.1145 

1 .3296E-0I 

1 -2370E-02 

-1 .593E-C4 

4.62728-02 

4.4767E+00 

0.2770 

8.00268*00 

-7.5276E  *00 

1 .1520E-OI 

0.9128 

0.3703 

8.8064 

1 .46008-01 

1 .33138-02 

-l .764E-04 

4.6504E-02 

1.35428*02 

0.2550 

8.8430E*S0 

-6. 1s86E*00 

1.07398-01 

0.9015 

0.3622 

7.1732 

1 .41108-01 

1 .47688-02 

-1 .225E-04 

4.7042E-02 

-2.84028+02 

0 . 2548 

8.0333E+OO 

-4.B1 79E+D0 

1 .DD02E-01 

0.8906 

0.3546 

6.0297 

1 .3650E-01 

1 . 59995-02 

-7.629E-05 

4. 55688-02 

-2.4612E+02 

0 . 2546 

8 . 9284E+00 

-2.56538+00 

0.9353E-02 

0.8769 

0 . 3436 

4.7932 

1 .2978E-0I 

1 .75098-02 

-2.237E-05 

4.3311 8-02 

-2 .5605E+02 

0.2543 

8.97358+00 

3. 46208-01 

7.95778-02 

0.0659 

0.3337 

3.9416 

1 .2405E-01 

1 .86878-02 

4.191 E-07 

4.13468-02 

-2 . 1 362E+02 

0 . 2536 

9.02428+09 

2.3924E+00 

7.47968-02 

0.0547 

0.3281 

3.4605 

1 .2219E-0I 

1 .9683E-02 

2.O40E-O5 

4.05428-02 

-1 .4J32E+02 

0.2533 

9.  07508*09 

4.05738+00 

7. 03758-02 

0.0503 

0.3230 

3.1020 

1 .2071E-01 

2.041EE-02 

5.987E-05 

4.0393E-02 

-9.0622E+O1 

0.2526 

9. 1 325E+09 

1 .6597E  *01 

5.85178-02 

0.0509 

D.3104 

1 .5452 

1 .6541E-01 

2.99218-02 

1 .OOOE-03 

5 . 9773E-02 

5.6225E+02 

0.2519 

9.16128*00 

1 .7335E+01 

5.60728-02 

0.0560 

0.3076 

1.5038 

1 .6462E-0I 

2.91 34E-02 

1.07IE-03 

5.90948-02 

2.5016E*02 

0.2515 

9.19008*09 

1 .31038*01 

5.36428-02 

0.0659 

0 .3048 

1 .4634 

1 .63S5E-01 

2.8295E-02 

1 .145E-03 

5. 83978-02 

2.7694E*02 

0.2512 

9.25578*00 

1 .9730E*  01 

4.09208-02 

0.6838 

0.2994 

1 .3849 

1.62648-01 

2.6524E-02 

1 .301E-03 

5.7092E-02 

2.8957E+02 

0.2505 

9.32148*00 

2 , 1 2408  +  01 

4.50008-02 

0.9018 

0.2948 

1  .3195 

1 .62008-01 

2.4896E-02 

1  .442E-03 

5. 60348-02 

2.8777E+02 

0.2498 

9.3972E*00 

2 .24I5E+01 

4.23298-02 

0.9161 

0.2914 

1 .2714 

1 .6237E-01 

2.3725E-02 

1 .552E-03 

5.55458-02 

1 .95648+02 

0.2491 

9. 4731 E*00 

2 . 351 9E  *  01 

4 . 0002E-C2 

0.9302 

0.2884 

1 .2289 

1 .6298E-01 

2.2636E-02 

1  .6518-03 

5.51658-02 

1 .9301E+02 

0.2464 

9.56158*00 

2.4370E+01 

3-84218-02 

0.9412 

0.2661 

1 . 1 964 

1 .64248-01 

2 . 1 075E-O2 

1 .727E-03 

5.5153E-02 

1 .2685E*02 

0.2477 

9. 6500E*00 

2.5103E+OI 

3.69708-02 

0.9522 

0.2039 

1.1664 

1 .65588-01 

2. 1 144E-02 

1 .7988-01 

5.5174E-02 

1  .26B3E  +  02 

0.2469 

9.75458*00 

2 .5659E+01 

3.6345E-02 

0.9562 

0.2824 

1.1476 

1 .6748E-01 

2.06348-02 

1  .8408-03 

5.5590E-02 

5.80108+01 

0.2460 

9. 8591 E*00 

2.61168*01 

3.57458-02 

0 . 9642 

0.2609 

1.1296 

1 .6938E-0I 

2 . 0528E-02 

1 .8608-03 

5 .60018-02 

5.8299E+01 

0.2451 

9. 98452*00 

2.6637E+01 

3.50658-02 

0.9714 

0.2792 

1.1095 

1. 71668-01 

2.01 77E-02 

1  . 923E-03 

5.64978-02 

5.736IE+01 

0.2441 

(i)  Best  fit  boundary  layer  solution. 
Figure  14.-  Continued. 


AEDCTH-BO-19 


ITERATION  NO.  13  (AVERAGE  OF  INVISCXD  AND  VISCOUS  SOLUTIONS  I 
PLOT  OF  CP  AT  EQUAL  XI-INCRENENTS 


I 

XB 

YB 

FN 

OS  CP 

1 

0.0 

0-0 

0.0 

0.0  1.1704  ♦ 

« 

2 

0.068 

0.017 

0.538 

0.695  0.5617 

♦ 

* 

3 

0.137 

0.034 

0.565 

0.727  0.5081 

+ 

* 

A 

0.207 

0.052 

0.642 

0.819  0.3472 

♦ 

* 

5 

0.2S0 

0.070 

0.665 

0.846  0.2967 

4 

* 

6 

0.356 

0.069 

0.679 

0.863  0.2663 

4 

* 

7 

0.435 

0.109 

0.689 

0.874  0.2446 

* 

X 

S 

0.519 

0.131 

0.697 

0.883  0.2274 

4 

« 

A 

0.608 

0.153 

0.704 

0.891  0.2123 

♦ 

X 

10 

0.703 

0.177 

0.711 

0.699  0.1976 

♦ 

ft 

11 

0.604 

0.202 

0.717 

0.907  0.1632 

+ 

X 

12 

0.911 

0.229 

0.724 

0.915  0.1676 

♦ 

* 

13 

1.026 

0.258 

0.732 

0.924  0.1503 

4 

* 

14 

t .  149 

0.289 

0.741 

0.934  0.1301 

4 

* 

15 

1.280 

0.322 

0.753 

0.947  0.1051 

4 

* 

16 

1  .418 

0.357 

0.767 

0.964  0.0720 

4 

« 

17 

1 .565 

0.393 

0.790 

0.969  0.0228 

4 

* 

IS 

1.720 

0.432 

0.840 

1.044  -0.0693 

4 

X 

19 

1  .685 

0.471 

0.956 

1.169  -0.3454 

♦  * 

20 

2.060 

0.496 

1.004 
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Figure  14.-  Continued. 
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Figure  14.-  Continued. 
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Figure  16.-  Two-dimensional  configuration  for  boundary- layer  calculation 
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Figure  18.-  Output  from  two-dimensional  boundary- layer  calculation. 
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NASA  DATA  COMPARISON  -  CONFIGURATION  1 
JET  SIMULATED  WITH  A  SOLID  STING 
BOATTAIL  L/D  =  0.8 
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B.  25 

0,477956 

(a)  First  50  cards. 

Figure  19.-  Input  data  for  inviscid  calculation. 
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(b! 

tfext  50  cards 

Figure  19. -  Continued. 
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(c)  Remaining  7  cards. 


Figure  19.-  Concluded. 
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NORMAL  COORD,  STRETCH  FOR  ALF=  1.500  - 
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Figure  23 


(a)  Normal  coordinates. 

-  Output  for  inviscid  calculation. 
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(b)  Tangential  coordinates. 


Figure  20.-  Continued. 
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0. 791 4E  *01 

0  7?3''E*01 

0  40Q4E*00 

0  0 

-0  871 Ct-01 

0.?’K6E-O1 

0 . 1098E*00 

49 

0  8051E+01 

0 . 7991 E*0 1 

0 .5000E*  00 

0.0 

-0.2599E*00 

0 . 3260E-0 1 

0.1 237E*00 

50 

0.6I46EMH 

0.6386E*01 

0 .497"*E  *00 

0.0 

-0.3447E+01 

0 . 7092E+00 

0 . 1409E*00 

51 

O.GE29E*01 

0 . 61 69E*  01 

0 . 490OE  *00 

0.0 

-0.6787E*01 

0. 7148E400 

0. 1618EX00 

52 

0.830 1 E  +  01 

0 .8241 E*01 

0.4796E+00 

0.0 

-0 . 9704E*0t 

0.73C9E400 

0.1866E*00 

53 

0.8363E*OI 

0 . 8303E  p  01 

0.4674E*00 

0.0 

“0. I22FE*02 

0. 7539E*00 

0 . 2 1 45F+00 

54 

0.64II»E*01 

0 . 83520*01 

0 , 4544E400 

0.0 

-0. 1451E*02 

0 . 7702E*00 

0.2626E+00 

55 

0 .8467E  *01 

O.64O7E*01 

0.440SE*00 

0.0 

-0. 1656E*02 

0.7926E+00 

0, 2647E+00 

56 

0.851 3E 1 01 

0. £4530*01 

0.4261 E  *00 

0.0 

-D.1850E*02 

0.8233E*00 

0. 2727E+00 

57 

0.85608*01 

0.£50Cr*01 

0.4097F*00 

a.o 

-0.2049F+02 

0.8550E+00 

0-P609E400 

53 

0.8610E*01 

0 ■ 8550E *0 1 

0. 389 9E*  00 

0.0 

-0.2264E+0E 

0 .  B942E+  l?0 

0.2396E+00 

59 

0  £6645*01 

0.86D5E*0I 

0.3658E*OQ 

0.0 

-0.2597E+0? 

0 , 9434E  +  DO 

0.2192E*00 

69 

0.87E4E*01 

0 . 8664E*0 1 

0.3362E*00 

0.0 

-  0 . 2771 E*02 

0.1053E*01 

0.1997E*00 

61 

0.67*?OE*01 

0.873CE*OI 

D  299-*E*0O 

0.  0 

-0.31 26E t02 

0.1I77E+01 

0. 1 81 2E*00 

$2 

0 . ES6 JE*01 

0 . 850  3E  *  0 1 

0 .2547E  *00 

0.0 

-0.2635E*OI 

“0.44"2E+02 

0 . 1 635E*00 

63 

Q.8Q'*3E*01 

0.8533E*01 

0  Z55QE+C0 

0.0 

-0.26G3E-04 

-0. I0I6E-0: 

0 . 1467E400 

64 

0 .9r3iE*01 

0 .8974E  *01 

0  2550E*00 

0.0 

-a.iaoaE-08 

0.2323E  06 

0. 1309E+00 

65 

0.91 J5E*01 

0.907EE+GI 

0.2550E*00 

0.0 

-0.0403F-09 

-0. 193JE-03 

0.1 I59E400 

66 

0 , 9250C  *01 

0.9199E+0I 

O.255OE+00 

0.0 

C.I146E-10 

-0.4979E- 1  1 

0.10I9E+00 

67 

0. 9731C >01 

0 , 9321 E  *01 

O.2550EIO0 

0.0 

-0.2712E-11 

-0.3707E-11 

O.S877E-0I 

66 

0.95J3EI01 

0.94730*01 

0 .2550E  *00 

0.0 

-0 . 7623E ■ 1 2 

o.i6ior-i2 

0.7654E-Q1 

69 

0 . 971 0E*  01 

0.9550E+01 

0.2550E+00 

0.0 

0.26C8E-12 

-0. 1067E-I3 

0.6522E-OI 

73 

0 . 991 9E*01 

0.9S5-E401 

0.25D0E*00 

0.0 

-0.  7049E-13 

-0.4717E- 1 5 

0.5490E-01 

71 

0. 101 7E*02 

0.1011 E*02 

0 . 2550C*00 

0.0 

0.1153E-13 

0.39I7E-I4 

0.4529C-0I 

72 

0. 1 048E  *02 

0. 1042E*02 

0 . 25 FOE*  00 

0.0 

0.9494E-15 

-0.32I5E-I5 

0.3669F-0I 

73 

o. ira6E*02 

o.iceoE+os 

0 . 2550E*  00 

0.0 

-0. 3641E-I 5 

-0.5504E-16 

0.2894E-01 

7<l 

0  11357*02 

0. 1 129E*02 

0.2550E»00 

0.0 

0. 2995E-1 6 

0.1040E-I6 

0 . 221 9E“01 

75 

0. 1 201 E*  02 

0.1I95E402 

0.2550E400 

0.0 

0.4933E-I7 

0 . 1  323E-17 

9.1630E-01 

76 

0 . 1 29JE  >92 

(J.12S7C*02 

0.£550E*00 

0.0 

O.D 

0.0 

B.1J32E-01 

77 

C. 1431E402 

0.1425E*02 

0,2550E*00 

0.0 

0.0 

0.0 

0.72*7E-0£ 

78 

0.16615*02 

0 . 165bE*02 

0 . 2550E  *0  0 

0.0 

0.0 

0.0 

0.4076E-02 

79 

0 . 21 21E  *02 

0.21 15E*02 

0.2550E*00 

0.0 

0.0 

0.0 

0.1812E-02 

80 

0  T50IE*02 

0 . 34 95E*  02 

0 . 2550E*0D 

0.0 

0.0 

0.0 

0.4530E-03 

81 

0.10000*31 

0. 1000E*31 

0.2550E400 

0.0 

0.0 

0.0 

0 . 1 0C0E-29 

(b) 

Concluded. 

Figure  20.-  Continued 
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ITERA’IOtl  HO.  1 


VO 

•jj 


If 

□  PriAX 

ID 

JD 

HHiX 

IR 

1 

0.405E-01 

61 

31 

0 . 31 2E<  02 

1 

8 

0 . ?37E  01 

2 

31 

0 . 1 19E402 

1 

3 

o.ieiE-oi 

1 

31 

0.I6£E*02 

1 

4 

0 . 1 64E-C1 

2 

31 

0 .6668*01 

59 

5 

0.I33E-01 

1 

31 

0 .1252*02 

1 

6 

0.967E-02 

58 

31 

0.446E+01 

58 

7 

O.845E-02 

S3 

31 

0  4ME.01 

62 

6 

0 , 757E-02 

57 

31 

0.3£4E*C1 

62 

9 

0  68 1 E - 02 

57 

31 

0.362E*01 

62 

ID 

0.61 7E-02 

56 

31 

0.366EI01 

1 

It 

0.562E-O2 

56 

31 

0. 3762*01 

62 

1  £ 

0.507E-02 

55 

31 

0 . 311 E  *01 

62 

13 

0.474E-02 

61 

31 

0,2981*01 

62 

14 

0. 454E-02 

61 

31 

0 . £862+01 

62 

15 

0.435E-02 

61 

31 

0.2/StiDI 

62 

16 

0.41 7E-02 

61 

31 

O.EfcSEtOI 

62 

17 

0.40  IE-02 

61 

31 

0.  ESSE  *01 

62 

18 

0.337E-02 

61 

31 

0. 296E  *01 

62 

19 

0. 373E-02 

61 

31 

0.2382*01 

62 

20 

0.36  IE-02 

61 

31 

O.230C+01 

62 

RHAX=  0.2JE*01,  CC3VR  =  0.15E-01 


JR 

I  SUB 

I  SUP 

RAVG 

FF1 

<?F3 

31 

1 

0 

0.149E+0O 

1.400 

0.100 

31 

1 

0 

0.1042*00 

1.400 

0.100 

31 

1 

0 

0.963E-01 

1.400 

0.1C0 

31 

I 

0 

0.B25E  01 

1  .400 

0.100 

31 

1 

0 

0.  788E  01 

1  .400 

0.100 

31 

1 

0 

0.70SE-0I 

1  .400 

0.100 

31 

1 

0 

O.637E-01 

1  .400 

0.100 

31 

1 

0 

0.603E-0! 

1  .400 

o.too 

31 

1 

0 

0.530E  01 

1  .400 

0.100 

31 

1 

0 

0.567E-OI 

1.400 

0.100 

31 

1 

0 

0.536E-0I 

1  .400 

0.100 

31 

1 

0 

0.5D7E-OI 

1  .400 

0.100 

31 

1 

0 

0.49IE-01 

1  .400 

0.  ICO 

31 

1 

0 

0.486E-01 

1  .400 

0.100 

31 

1 

0 

0.476E-01 

1  .400 

0.100 

31 

1 

0 

0.46IE-O1 

1  .400 

0.100 

31 

1 

0 

0.447E-01 

1  .400 

0 . 1  DO 

31 

1 

0 

0.435E-01 

1.400 

0.100 

31 

1 

0 

0.431E-0I 

1.400 

0.100 

31 

1 

0 

0.423E-0I 

1  .400 

o.ioo 

MS 


0 

0 

a 

o 

o 

o 

o 

0 

0 

0 

0 

a 

o 

o 

o 

s 

4 

4 

5 
5 


(C) 


Iteration  history 


SEC/C1C 

a.  iso 
o.  iso 
a.  ion 
0.180 
o.teo 
0.160 
0.160 
0.180 
0 . 1 60 
0.180 
0.160 
0.180 
0.180 
0.180 
0.180 
0.180 
0.170 
0.170 
0-180 
0.180 


Figure  20.-  Continued 
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ITERATION  MO.  1 


PLOT  OF  CP  AT  EQUAL  XI-INCREMENTS 


I 

XB 

YB 

Ffl 

QS  CP 

1 

0.0 

0.0 

0.0 

0.0  1.1704  + 

ft 

z 

0.068 

0.017 

0.835 

1 .038  -0.0773 

♦ 

ft 

3 

0.137 

0.034 

0.590 

0.758  0.4557 

+ 

ft 

4 

0.207 

0.052 

0.641 

0.818  0.3479 

♦ 

* 

5 

0.280 

0.070 

0.663 

O.B44  0.3004 

+ 

ft 

6 

0.356 

0.039 

0.676 

0-859  0.2733 

4 

ft 

7 

0.435 

0.109 

0.694 

0.869  0.2548 

♦ 

ft 

S 

0.519 

0.130 

0.691 

0.877  0.2404 

4 

ft 

9 

0.608 

0.152 

0.697 

0.883  0.2279 

4 

w 

10 

0.703 

0.176 

0.702 

0.890  0.2159 

♦ 

* 

11 

0.804 

0.201 

0.708 

0.896  0.2033 

♦ 

K 

12 

0.912 

0.228 

0.714 

0.933  0.1697 

♦ 

ft 

13 

1.027 

0.257 

0 . 721 

0  911  0.1741 

4 

ft 

14 

1.149 

0.287 

0.730 

0.921  0.1554 

♦ 

ft 

15 

1  .280 

0.3SO 

0-740 

0.933  0.1320 

4 

ft 

16 

1.419 

0.355 

0.755 

0.949  0.1006 

* 

ft 

17 

1.566 

0.392 

0.776 

0.973  0.0532 

♦ 

ft 

18 

1.721 

0.430 

0.823 

1.025  -0.0506 

4 

ft 

19 

1 .885 

0.470 

0.943 

1.154  -0.3144 

4  ft 

20 

2.060 

0.495 

1  .003 

1.215  -0.4407 

4 

21 

2.246 

0  500 

0.903 

1.117  -0.2375 

4 

ft 

22 

2.440 

0.500 

0.841 

1.045  -0.0903 

4 

ft 

23 

2.642 

0.500 

0  825 

1.027  -0.0550 

* 

ft 

24 

2  851 

0.500 

0  816 

1.016  -0. 0260 

♦ 

ft 

25 

1.088 

0.5C0 

O  811 

1.012  -0.0245 

4 

ft 

26 

3.290 

0.500 

0.808 

1.009  -0.0172 

4 

ft 

27 

3.518 

0.500 

0.806 

1.006  -0.0123 

4 

ft 

28 

3.751 

0.500 

O.B04 

1.004  -0. 0090 

♦ 

ft 

29 

3.988 

0.500 

0.803 

1.003  -0.0066 

♦ 

ft 

30 

4.229 

0.500 

0.002 

1.002  -0.0049 

4 

ft 

31 

4.471 

0.500 

0.802 

1.002  -0.0037 

4 

ft 

IS 

4.715 

0.500 

0.801 

1.001  -0.0023 

4 

ft 

33 

4.958 

0.500 

0.801 

1.001  -0.0321 

4 

ft 

34 

5.201 

0.500 

0.B01 

1 .001  -0.0016 

♦ 

ft 

35 

5.442 

0.500 

O.SOI 

1 .001  -0.001 1 

4 

ft 

36 

5.679 

0.500 

0.800 

1 .000  -0.0309 

4 

ft 

37 

5.912 

0.500 

0.800 

1 .000  -0.0007 

4 

ft 

38 

6.139 

0.500 

0.800 

1  .ODD  -0.0098 

4 

ft 

(d)  Tabulated 

solution  and  plot 

of 

Figure  20.-  Continued. 
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* 


o 

is* 


39 

6.359 

0.500 

0.801 

1.001  -0.0012 

90 

6.572 

0.500 

0.801 

1.091  -0.0018 

91 

6.776 

0.500 

0.801 

1.001  -0-0026 

92 

6.970 

0.500 

0.802 

1  .002  -0  .0036 

93 

7.153 

0.500 

0.802 

1.002  -0  0049 

99 

7.324 

0.500 

0.804 

1.004  -0.0083 

95 

7.483 

0.500 

0.603 

1.008  >0.0163 

96 

7.630 

0.500 

0.815 

1.017  -0.0340 

97 

7.763 

0.500 

0.82? 

1.031  -0.0628 

98 

7.584 

0.500 

0  5=1 

1.C55  -0.1113 

99 

7.991 

0.509 

0.  912 

1.122  -0.2473 

50 

8.066 

0.497 

0.973 

1.202  -0.4190 

51 

8.  169 

0.490 

1.031 

1  .234  -0.4935 

O.P'.I 

0.430 

1  .040 

1 .234  -0.5170 

53 

8.303 

0.467 

1.033 

1.217  -0.5037 

54 

6. 358 

0.454 

1.013 

1.186  -0.4623 

55 

8.407 

0.441 

0.989 

1.151  -0.4111 

56 

8.453 

0.426 

0.963 

1 .114  -0.3565 

57 

8.500 

0.410 

0.973 

1.061  -0.2712 

58 

8.550 

0.390 

0,868 

0,992  -0.1509 

59 

8.605 

0.366 

0.794 

0.900  0.0134 

60 

8,664 

0.336 

0.685 

0.770  0.253! 

61 

6.730 

0.299 

0.490 

0. 544  0.6544 

62 

8.803 

0.255 

0.467 

0.632  0  6601 

63 

8. 833 

0.255 

0.561 

0.722  0  5153 

64 

8.974 

0.255 

0.627 

0.801  0.3783 

65 

9. 075 

0.255 

0.673 

0.856  0.2791 

66 

9.190 

0.255 

0.703 

0.8=6  0.2029 

67 

9.321 

0.255 

0.735 

0 . 92  7  0.1444 

68 

9.473 

0.255 

0.755 

0,949  0.1002 

69 

9.650 

0  255 

0.770 

0.966  0.0676 

70 

9.859 

0.255 

0.780 

0.978  0.0441 

71 

10.110 

0.255 

0.787 

0.536  0.0277 

72 

10.417 

0.255 

0.772 

0.992  0.0167 

73 

10.809 

0.255 

0.796 

0.995  0.0096 

74 

11.293 

0.255 

0.743 

0.997  0.0051 

75 

11.950 

0.255 

0.799 

0.999  0.0025 

76 

12.870 

0.255 

0.799 

0.999  0.C011 

77 

14.250 

0.255 

0.800 

1.000  0,0004 

78 

16.550 

0.255 

0.600 

1.000  0.0001 

79 

21 .150 

0.255 

0 . 600 

1.000  0.0100 

60 

34.943 

0.255 

O.BOO 

1.000  O.0 

81. 

m****** 

0.255 

0.090 

1,000  0.0 

TOTAL  BOTY  DRAG  COEFFICIENT*  Oi 07352 


AFTERBODY  DRAG  COEFFICIENT*  0.04209 


(d)  Concluded. 
Figure  20.-  Continued 
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AEDC-TR-80-19 


ITERATION  NO.  I 


WO.  OF  SONIC  RTS.  OH 

EACH 

CONS ANT -J 

LINE, 

STARTING  WITH 

JtJMAX 

(BOOT  SURFACE* 

AT 

LEFT 

4  0  0  0  0  0  0 

0  0 

0  0  0 

0  0 

0  0  0  0  0 

0  0 

D  0  0  0  0 

0 

0  0 

H= 


4 


XS(K  ),  K-* , >H 

0.205E+01  0.2C7E*D1  0.610EK)t  0.638E+O1 


TStKI.  K-l . N 

0.495E«00  0 .496  Eh  CO  0.49iE*0D  0.447E*0D 


(e)  Sonic  point  distribution. 


Figure  20.-  Concluded 
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TABLE  I.  RELATION  BETWEEN  EXTERNAL  DATA  SETS 
AND  INPUT/OUTPUT  LOGICAL  UNIT  NUMBERS 


Logical  Unit 

Data  File 

(fig.  12)  Input  Output 


1  14  9 

2  13  8 

3  12  2 

4  15  10 

5  11  3 


97 
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TABLE  II.  BASIC  INPUT  DATA 

Item 

No.  Variables 

1  TITLE  (three  cards) 

2  NRSTRT, N3, ILIM 

If  NRSTRT  /  0  and  N3  f  0,  skip  to  item  12 

If  NRSTRT  /  0  and  N3  /  0  stop.  No  further 
input  is  necessary 

3  LPROG, N1,N2, IBL, IUNIT, MIT 

4  GAM, AMINF 

If  LPROG  =  1  skip  to  item  13 

5  IXY 

If  IXY  =  0r  skip  next  card 

6  XO,YO  (IXY  cards) 

7  IMAX, JMAX , MHALF , KLOSE , LREADP 

8  DNDZO , XIXM, XM, DSDXIM , XBT , DMAX , XZNEW 

9  PLUMIN 

If  PLUMIN  =  FALSE,  skip  next  card 

10  GAMAP , PTPPFS ,  AMP , THETAP , TTP , GMP , X J , GC 
If  LPROG  =  -1,  skip  remaining  cards 

11  LSEP, LITER, LNSHK 

If  LSEP  =  FALSE,  skip  next  card 

12  XSEP , DTHET 

If  NRSTRT  ^  0  stop.  No  further  input  is 
necessary 


Format 

20A4 

315 


615 

15 

2F10. 0 
515 

7F10.0 

L5 

7F10.0 

L5 

2F10 . 0 
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TABLE  II.  CONCLUDED 

Item 


NO. 

Variables 

Format 

13 

IOPT ,  K , LVAR1 , LSHAPE ,  LIC  , LDSTAR , LSHPBL 

If  LPROG  ^  1  or  if  LPROG  =  1  and  both 

IOPT  =  1  and  LVARl  =  2,  skip  items  14 
and  15 

715 

14 

NVAR 

15 

15 

XVAR , VAR 

2F10.0 

16 

EL , PT , TT , TWONTT , VI SC , RGAS , SCON , DFACT 

8F10.0 

17 

XZ , RLEN , XT , DXP , HLIM 

If  LSHAPE  ^  0,  skip  next  card 

If  LIC  =  1,  input  CFCl  and  DELTAl 

If  LIC  =  2 ,  input  CFCl  and  DELSTl 

5F10.0 

18 

CFCl, DELTAl  {or  DELSTl) 

If  IOPT  =  1,  skip  next  card 

2F10 . 0 

19 

UEl , DUEDX 

If  LSHPBL  =  0  and  LPROG  =  0,  skip  items  20 
and  21 

2F10.0 

20 

NR 

15 

21 

XRP , RL 

2F10.0 

99 
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TABLE  III.  SUMMARY  OF  INPUT  DATA  FOR  RESTARTING 

Item 

No .  Variables 

1  TITLE  (Three  cards) 

2  NRSTRT , N3 , ILIM 

If  N3  =  0  stop.  No  further  cards  are 
necessary 

3  XSEP , THETS 


Format 

20A4 

515 

2F1Q.0 
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H 


i 


M 


m 


P 

r 

S 


s 

T 


u,  v 


x,r 


6 


Y 


<5 


* 

i 


0 


s 


LIST  OF  SYMBOLS 

coefficients  in  eq.  {24) 
speed  of  sound 
term  in  eq.  {24) 
width  of  mixing  layer  {fig.  2) 
skin  fiction  coefficient 
maximum  body  diameter 
quantity  defined  in  eq.  (14) 
transfomred  shape  factor,  5* / 6 ^ 

Mach  number 

Ill  m2 

2 

pressure 

radius 

VTte  '  1 

root-mean- square  error,  eq.  (25) 
temperature 

x,r  velocity  components 

physical  coordinates  defined  in  figure  2 

location  of  separation  point  (fig.  1) 

eddy  viscosity  factor,  T/y(5u/3y) 

ratio  of  specific  heats 

transformed  boundary- layer  thickness 

displacement  thickness  5 

transformed  displacement  thickness,  J  (1  -  u/ue)dy 

o 

half-angle  of  conical  displacement  surface 
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6 


i 


u 


v 

p 

T 

t 

to 


LIST  OF  SYMBOLS  (Concluded) 

6fi  u 

transformed  momentum  thickness,  J  — 
molecular  viscoity 
u/p 


o  ue 


density 

shear  stress 

angle  in  eqs.  (4 ) — ( 7 ) 

underrelaxation  factor  [eq.  (2)] 


dy 


Subscripts 

e  refers  to  boundary- layer-edge 

f  refers  to  conditions  after  a  shock  wave 

i  refers  to  conditions  preceding  a  shock  wave,  or  con¬ 

ditions  at  the  inner  edge  of  an  exhaust  jet  mixing 
layer 

I  refers  to  inviscid  flow 

j  refers  to  conditions  in  exhaust  nozzle 

o  refers  to  the  undisturbed  free  stream 

s  refers  to  separation  point 

t  denotes  stagnation  point 

V  refers  to  viscous  flow 

w  refers  to  conditions  on  a  solid  surface 


Special  notation 

(  )  refers  to  differentiation  with  respect  to  x 

X 
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